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, Abstract 

O^. Lattice quantum field theory is a well established branch of modern quantum 

field theory (QFT). However, it has only peripherally been used for the in- 
vestigation of Casimir systems - i.e. for systems in which quantum fields are 
distorted by their interaction with classical background objects. This arti- 
5^ I cle presents a Hamiltonian lattice formulation of static Casimir systems at a 

level of generality appropriate for an introductory investigation. Background 
structure - represented by a lattice potential V{x) - is introduced along one 
spatial direction with translation invariance in all other spatial directions. It 
is simple to extend this formulation to include arbitrary background structure 
in more than one spatial direction. Following some general analysis two spe- 
cific finite ID lattice QFT systems are analyzed in detail. The first has three 
Dirichlet boundaries at the lattice sites x = 0, I and L (L > I > 0) with 
vanishing lattice potential V{x) everywhere in between. The vacuum energy 
and vacuum stress tensor T^'^ for this system are calculated in < rr < L. 
Very careful attention must be and is given to renormalization in the (con- 
tinuum) limit of vanishing lattice constant. Globally and locally this lattice 
system is seen to closely mimic the corresponding ID continuum system - as 
one would hope. Then we introduce a lattice potential V{x) = c/{x — xq)^ 



centered at x = xq < to the left of the boundary at x = and extending 
through this boundary and the middle Dirichlet boundary at a; = Z out to the 
right-hand boundary x = L > I and beyond. The vacuum energy and T^'^ 
arc calculated for this far more complicated system in the region < .x < L, 
again with very good results. The internal consistency of the lattice version 
of this system is carefully examined. Our conclusion is that finite-lattice for- 
mulation provides a powerful and effective tool, capable of solving completely 
many Casimir systems which could not possibly be handled using continuum 
methods. This is precisely our reason for introducing it. Future investigations 
(in one and more dimensions and in dynamical as well as static contexts) will 
display more fully the power of this method. 
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1 Introduction 



Let us here define Casimir QFT (quantum field theory) to be the study of quantum 
fields in flat space, fields which coexist with and interact with classical background 
objects occupying the same space (in the simplest case Dirichlet or Neumann bound- 
aries). This interaction distorts the quantum fields and all local observables con- 
nected with them away from the spatial uniformity (or translation invariance) that 
would characterize these fields in empty fiat space. This nonuniformity in fiat space 
caused by background objects is generally very pronounced, becoming even infinite 
as Dirichlet or Neumann boundary surfaces are approached. The background ob- 
jects themselves experience measurable back forces - a subtle example of Newton's 
third law. These back forces are called Casimir forces and their presence is often 
identified with the Casimir effect. However the deeper essence of the Casimir effect 
is the distortion of the quantum fields caused by the background objects. Viewed 
this way the Casimir effect is obviously an extremely fundamental and general phe- 
nomenon in quantum field theory. Any nonuniform classical background whatsoever 
will produce a Casimir effect. 

To the authors' knowledge all published work in Casimir QFT employs continuum 
QFT. In certain other branches of QFT, of course, lattice formulations have be- 
come commonplace. Quite well-developed methods exist for doing QFT on a spatial 
or spacetime lattice (see e.g. the book by Rothe |]1[). Generally in lattice-QFT 
applications spatial uniformity is assumed. Either one chooses the spatial lattice 
to be infinite, or one imposes periodicity on a finite spatial lattice. There are no 
background objects on the lattice. Rather, one is using the lattice to learn about 
nonperturbative properties of interacting fields which are difficult if not impossible 
to investigate in continuum theory. 

Some years ago two of the authors undertook the reformulation of Casimir QFT in 
lattice language. From the outset we have had two main objectives. 
First, to gain a sense of how well lattice QFT works when applied to Casimir prob- 
lems. By investigating a number of systems whose continuum versions can be solved 
explicitly we found that lattice QFT is able to mimic and reproduce in a quite de- 
tailed way both global and local features of continuum Casimir systems. While not 
unexpected this did need to be established quantitatively before we could proceed 
to our second and more important objective. 

That objective is, to use lattice QFT to solve Casimir problems which cannot be 
dealt with (adequately, or even at all) using continuum QFT. It is unfortunate that 
practically all Casimir systems of physical interest belong in this latter difficult-to- 
access category. 

On a lattice one can work both analytically and numerically. Analytic methods tend 
to rather closely parallel continuum theory and, while interesting, do not necessar- 
ily offer greatly increased computational power, especially on an infinite lattice. 
Numerical methods used on a finite lattice can be, however, very powerful. One 
can introduce static background objects on the finite spatial lattice and make the 
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quantum field interact with these objects, the interaction being described by 
a potential in a functional Schrodinger equation, whose stationary solutions are ex- 
pressible in terms of the eigenfunctions and the spectrum of an ordinary Schrodinger 
equation associated with the functional one. These spatial modes and the spectrum 
can be found numerically. Given these basic ingredients one can then compute ev- 
erything numerically - the vacuum energy and the stress-energy-momentum tensor, 
Casimir forces and so on. This can surely be done for a vast array of systems which 
could not be solved in continuum theory. 

Of course renormalization must be performed to make lattice observables physically 
meaningful. Everything one calculates on a finite lattice is finite. However, lat- 
tice quantum variables like vacuum energy and the vacuum expectation value of 
the stress-energy-momentum tensor T'^^ contain the seeds of ultraviolet (UV) di- 
vergences, i.e. terms which grow without limit in the continuum limit a — » (a 
is the lattice constant or spacing between lattice points). For finite and decreasing 
a these essentially informationless terms eventually overwhelm the physical signal: 
their systematic removal is utterly essential. Lattice quantities which vanish as 
a —>■ take care of themselves, of course. What remains is the "physical signal" - 
a contribution independent of a which therefore neither blows up nor vanishes as 
a ^ 0. All of this is the process of renormalization. 

For a number of years we have pursued this numerical strategy, for both standard 
and unconventional Casimir systems. To date none of this lattice QFT work has 
been reported (excepting the brief descriptions 0, 0]). We begin here the systematic 
development of lattice QFT with background objects and its application to nontriv- 
ial Casimir systems. Before describing the contents of the present article we pause 
to mention some nonstandard aspects of Casimir QFT. 

Traditional Casimir systems have "hard" boundary surfaces - i.e. boundary sur- 
faces dAi with precise spatial location on which the quantum field (f){x) has to 
satisfy some boundary condition. A perhaps more physically realistic boundary can 
be fashioned by attaching to a Dirichlet boundary dAi (on which (f){x) has to van- 
ish) a static potential V{x) which grows to infinity as x approaches dAi and falls 
off to zero more or less rapidly away from dAi. The interaction of this classical 
potential with the quantum field 0(x) is for simplicity assumed to be described by 
an interaction term quadratic in 0(x) : (1/2) 0^(x) V{x). We refer to such a modi- 
fication of a hard Dirichlet surface as semihardening 0, ^. The potential function 
V{x) represents (one can say) surface texture attached to dAi. The wave equa- 
tion [Oq — A + m^] = in usual hard-boundary scalar-field Casimir QFT with 
appropriate boundary conditions on is replaced by 

[(9^ - A + + V{x)] (j){x) = . (1.1) 

Here V{x) should be viewed as a (rather huge) generalization of ordinary bound- 
ary conditions. This potential can represent many kinds of background structure 
which interact with and are (for V{x) < oo) partially transparent to 0. The 
semihardening potentials above are just one kind of background spatial structure 
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among many which can be introduced into QFT (see e.g. refs. ^, |T^, |T5|). For us, 
Casimir QFT comprises all possible background structures interacting with quan- 
tum fields - a truly vast subject. Many workers, of course, use the term "Casimir 
effect" mainly in connection with boundary surfaces (see e.g. the reviews p!6|). 



Diffuse background structures in QFT are then often called "background fields". 

As mentioned earlier the fundamental phenomenon in Casimir QFT is the dis- 
tortion of otherwise spatially uniform fields and observable quantities away from 
spatial uniformity by the background V{x). To really understand a Casimir system 
one must work locally. There is no better physical quantity to compute than the 
vacuum expectation value of the stress-energy-momentum tensor (hereafter stress 
tensor) of a Casimir system which we will denote by {T^^{x)) in the following. 
This quantity reveals the spatial nonuniformity of the system, its nonuniform vac- 
uum energy density and its Casimir forces. {T^^{x) ) can be computed from the 
spatial modes and spectrum of 0(x), both in continuum theory and on a lattice. 
When working numerically on a finite lattice (using a lattice potential V{x)) one 
can obtain (T^'^(x) ) numerically for practically any kind of background structure. 



In sec. 2 we begin our analysis of Casimir systems of the following type: A space 
of dimension d is chosen with free (or periodic) "boundary conditions" in each of the 
directions X2, ■ ■ ■ , Xd- Along the xi axis we place parallel Dirichlet boundaries at 
Xi = and xi = I. In sec. 2 we do not concern ourselves with the regions xi < 
and Xi > I external to the interval < Xi < / of interest. These external regions 
can be "attached" later on and given suitable global structure (e.g. periodicity in 
Xi over some larger interval) to complete the specification of the physical system. 
We postpone doing this without loss of generality. In the interval < Xi < I 
we introduce an arbitrary background potential V{xi). The formulation of lattice 
scalar QFT on this spatial manifold is reviewed in detail in sec. 2, where detailed 
lattice formulae for (T^'^) are given and the ground is prepared for our very specific 
lattice analyses to follow in sees. 3, 4. 

In sec. 3 we go to one dimension ci = 1, set the background potential V{x) = 
and study scalar lattice QFT on a ID lattice < x < Z with Dirichlet walls 
at X = 0, / and no structure in between. Choosing lattice constant a means 
/ = Na , X = na where n = 0, 1, ■ ■ ■ , labels the lattice sites. The continuum 
limit a — >■ (with / fixed of course) means N = l/a ^ oo. As described above we 
must identify and subtract those terms in the lattice vacuum energy and {T^'^) of 
this system which diverge for a — 0. The remaining a-independent parts of these 
functions are then to be brought into comparison with continuum theory. In sec. 3 
we present detailed analysis of this continuum limit on the ID lattice. To enable 
us to make definite statements about the Casimir force on the boundary x = / we 
introduce an "external" region I < x < L with L S> / which can be analyzed on 
the lattice exactly as we do the "internal" region < x < /. We then have a more 
complete physical system than just the internal region < x < /. To make it a re- 
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ally complete physical system the region x < needs to be included and the global 
structure of the x axis must be specified. However there is no need to do all this 
here. Our goal is to show how well hard-boundary renormalized lattice QFT mimics 
continuum QFT. That goal is nicely achieved using the two adjacent ID lattices 
0<x<l, l<x<L. 

Sec. 4 tackles a more difficult problem. Keeping spatial dimension d — Iwe intro- 
duce the Bessel semihardening potential V{x) — cj {x — x^f' centered at x = < 0, 
allowing this potential to extend arbitrarily far into x > 0. Immersed in Vix) we 
position Dirichlet point boundaries at a; = 0, /, L as before. One can find the lattice 
vacuum energy E'""''^ and (T'*'^) numerically for this system with arbitrary back- 
ground potential Vix)- In sec. 4.1 we show how to renormalize E""^ and (T^*^) for 
arbitrary V{x). Then we adopt the Bessel potential above and begin our principal 
numerical analysis. In a sense the numerical work in sec. 4.2 is the central part 
of this article. Using E'"^.'^'^^ we calculate globally the Casimir force on the middle 
boundary x = I due to the other Dirichlet boundaries and the potential V{x). Us- 
ing {T'^^)ren we Calculate locally the same Casimir force. Very precise agreement is 
found. Some other things are computed as well, but let us limit our description of 
sec. 4 to the comment just made. This agreement between global and local calcu- 
lations of a physical Casimir force in a quite nontrivial background potential gives 
us great confidence in the methods we are using. We have no doubt that equally 
good results can be gotten for practically any background V[x). Indeed we have 
had comparable success with other potentials. 

Concluding remarks are made in sec. 5. There we indicate the directions in lat- 
tice QFT theory we presently are pursuing. Three appendices dealing with technical 
matters enable us to present a substantially condensed narrative. 



2 Scalar field on a finite lattice with background 

To construct our Casimir system we begin with fiat d-dimensional space and give 
it background structure along just one direction, say Xi. The other directions 
x±_ = (x2, • • • , Xd) are free but - in preparation for the finite lattice - we make 
each of these coordinates periodic in —L < Xj < L with period 2L. Periodicity is, 
of course, closely akin to "free boundary conditions" , neither of these being a true 
boundary condition at all. Thus x± specifies position on a spatial -dimensional 
torus T'^~^ representing the free directions and x = (xi, x±) is position in space. A 
real scalar field (j){x) with mass m is now defined on this space. 0(a;) interacts with 
and is made nonuniform by and exerts a back force on the background structure 
along the Xi direction. This interaction is what we wish to study in detail. 
For completeness let us mention something we do not wish to study more than nec- 



4 



essary: the so-called "topological Casimir effect". When a quantum field is 
defined on a topologically nontrivial space, is forced to be (nontrivially) dif- 
ferent than it would be in free infinite space. This is itself a kind of Casimir effect. 
The vacuum energy density of 0(x) is shifted, vacuum stresses or Casimir forces are 
exerted by the modified field directly on the spatial manifold causing the modifica- 
tion etc (see e.g. the reviews |TB[)- When, as here, one defines one's field on a 
(d-l)-dimensional torus T'^-^ (the periodic directions X2, ■ ■ ■ , Xd) there arises a topo- 
logical Casimir effect parametrized by the finite circumferences of T^~^. However, 
no "objects" have been introduced anywhere; periodic boundary conditions do not 
involve real boundaries. The T'^~^ topological Casimir effect smoothly disappears 
as the torus circumferences become infinite. 

Topological Casimir effects caused by global spatial structure coexist with the field 
distortion caused by locally nonuniform spatial backgrounds, or objects in space. 
One should, however, make the effort to keep these two aspects of QFT distinct 
in one's thinking. This seems particularly important when one replaces continuum 
QFT by lattice QFT. For numerical work the lattice has to be finite, inevitably 
introducing the topological effect under discussion, which lattice workers and others 
usually call "finite-size effects" . Unless one wishes specifically to study them, finite- 
size effects tend for obvious reasons to be regarded as a nuisance. 
Many readers may be familiar with finite-size scaling theory (see e.g. refs [0, |19[) 
in the general area of statistical systems. This and other related work one finite-size 
effects in QFT has rather limited overlap with and relevance for the kind of investi- 
gation we are concerned with here. Our focus is on the distorting effect background 
objects have on quantum fields. Finite-size scaling theory usually deals with bulk 
effects and large scale collective phenomena near critical points, and the effect of 
small system size on such phenomena. Phase transitions are involved, and tempera- 
ture plays an essential role. In Casimir QFT temperature is something of a sideshow 
(although not completely without interest), and phase transitions normally play no 
role whatsoever. True, the effects of surfaces and corners in statistical models are 
sometimes taken into account (see e.g. ^^), and this overlaps with our work more 
substantially. However, an attempt to pursue this would take us far outside the 
scope of the present work, and necessitate substantially increased length in an al- 
ready long paper. 

Resuming now our mathematical discussion, we position along xi two planar Dirich- 
let boundaries at Xi = and Xi = I. Between these boundaries we subject the 
quantum field (f){x) to a background potential V{xi) > 0. Lacking external regions 
Xi < and Xi > I this Casimir system is physically incomplete. However it is easy 
to attach these external regions. As long as one does not try to discuss Casimir 
forces on the boundaries xi = 0, / there is no need to be concerned about anything 
outside < X < I. Our purpose in this section will be the detailed formulation of 
the lattice quantum theory in < x < Z in the presence of arbitrary background 
potential V{xi). We begin with a brief summary of the continuum QFT we wish to 
put on a spatial lattice. 
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2.1 Continuum system 

The classical Lagrangian of the real scalar field interacting with a background, rep- 
resented by a classical potential V{xi) becoming infinite at xi = 0, /, is 

Jd'^xC ^ Jd'^x^ - (m^ + V{xi))(f>^) (2.1) 

= Jd'^x^ (02 - (j){m^ + V{xi) - A)4)) . 

Here in the partial integration used to reach the final equality the surface term 
vanishes because of the boundary conditions imposed on the classical field (jy^x^t): 

0(a;i = 0,f^,t) = 0(a;i = /,f^,t) = (due to ^(0) = V(/) = oo) (2.2) 
(j){xi,x^,t) — (l){xi,x± + 2Lei,t) , i — 2,...,d (periodicity), 

6*1 being the unit vector in i-direction. Using the canonically conjugate field 

o(p(x, t) 

we obtain from L by Legendre transformation the classical Hamiltonian 

H = J d'^x7r{x,t)<p{x,t) - L = J d'^x^(7r'^ + (j){m^ + V{xi) - A)(j)) . (2.4) 

Canonical quantization means replacing the classical field 0(x, t) and its conjugate 
momentum 7r{x, t) by hermitian operators $(x, t) and Il{x, t) which obey the equal- 
time commutation relations 

mx,t),^x',t)] = [Uix,t),U{x',t)] = 0, (2.5) 
mx,t),U{x',t)] = iS{x-x'). (2.6) 

Following standard procedure one can realize these commutation relations by ex- 
panding the field operator in terms of creation and annihilation operators a|. and d^ 
which satisfy [afc,a|./] = 6k^k' and [dk,dk'] = [a|.,aj.,] = ; 

4>{x,t) = Yl Idke-'^^'ukix) + dW^'^'ulix)] . (2.7) 

Here the Uk{x) and are eigenfunctions and the associated eigenvalues of the op- 
erator [—A + m'^ + V{xi)] , i.e. 

[m^ + V{xi) - A]u\x) = elu''{x) . (2.8) 
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Separating variables 

u''(x) ^v'''(xi)w''^{x2)---w''''{xd) = v'''(xi)w''^(xj_) , (2.9) 



k ^ {ki,k2, . . . ,kd) = {ki,k±) , 



leads to 



V(xO - v'^ix^) = p^^^^H^i) , (2.10) 

-7^-T^w'='(xi) = u;^.w'^'(xi) , i = 2,...,ci and (2.11) 
4 = rn' + pl+Y^ujl , (2.12) 



1=2 

where the boundary conditions on v''^ and w^' are 



v*^i(0) = v'^'il) = , (2.13) 
w'''{xi) = i(;''^(xi + 2L) , i^2,...,d. 



Obviously for the free directions 



with eigenvalues 



w''{xi)^J^e'i''^-^ , (2.14) 



< = ^kl heZ. (2.15) 



The functions v''^{xi) of course depend on the potential V{xi). These functions are 
chosen to be real. 
Altogether we have 



d-l 

u'^ix) = i;'^(xi)u>^^(fx) =t'''Ha;i) (^— ) ' e'^^^^^ , (2.16) 



d ^2 



Sk ^ m . , ^ 

j=2 ^ 



The set of modes {u''{x)} is complete and orthonormal 

Y,u\x)u^*{x') = 6{x-x') and j d'^xu^{x)u^'* {x) = 5k,k' , (2.18) 



where 



oo oo 



H ■■■ X] ••• > h,k' ^ Ski,k[Sk.2,k!,---hd,k'^ ■ (2.19) 

fe A;i=l fc2=— oo kd=—oo 
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Due to periodicity in X2, ■ ■ ■ , we need to introduce here the periodic 5-function 
S{x) along these directions, distinguished from the usual 5-function for a noncompact 
direction by its tilde. The spatial ^-function is then 

S(x-x') = S{xi-x[)5{x±-x'^) . (2.20) 

Our attention now shifts to the vacuum stress (energy-momentum) tensor of the 
system. The classical canonical stress tensor for a real scalar field is 

T^- = l((9»(9» + (9»(9»)-(?^^£ , (2.21) 

with £ given by the first equality in eq. (p.l|). Again standard procedure (see e.g. 
[IC]) leads to the following mode-sum formulae for the vacuum expectation value of 
the operator T^'^: 

(0|T»»(x)|0) = E '^"^Lr-i (2-22) 



, Ask dxi \ ' dxi 

i2el-pl + V{x^))v'^{x,)' 

k 



{0\T'\xm = E 



^ 1 /dv'^^^xiY 
+ ^ Ask (2L)^-i 1^ dx, ^ 

{Pl-V{x,))v^^{x,f 



2ek {2LY-' 



and for i = 2 d 



^ {pl-Vix,))v'^^ix,f 1 ( dv''^{x,) y 

^ 4ek{2LY-' ^^4eki2LY'^[ dx, ) 

{0\T\x)\0) = y ^fc. (2.24) 
^ ' ^ ^' ^ ^2£fc(2L)'^-i ^ ^ 

V 4£fc (2L)^-i dx, ^""'^ dx, j 
{pl-V{x,) + 2ujl)v^^{x,f 



E 

k 

E 



^eu{2LY-^ 
( dv^^{xiy 



k 



Ask (2L)^-i V -9x1 
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The first equality in eq. ( |2.22| ) shows that 
as one expects. 



(2.25) 



2.2 Lattice system 

The Casimir system will now be redefined on a finite lattice, which is given the 
same lattice constant a in all spatial directions. Thus spatial position becomes x = 
(xi, ■ • • , Xd) = a (^1, ■ ■ ■ , where rii, ■ ■ ■ , nd all take integer values. The planar 
Dirichlet boundaries at Xi = 0, / have lattice positions rii = and rii = l/a = Ni. 
The continuum potential V{xi) becomes a lattice potential V{ni) > 0. The free 
directions x± = an± are still periodic with period 2L: 



— L < Xj < L 



L 

<ni< 

a 



L 



AT, 



, d . 



(2.26) 



i.e. the lattice is periodic with period 2N_[_ in the spatial directions i = 2, ■ ■ ■ , d. 
On the lattice the Dirac ^-function 6{xi — x'l) is replaced by Sn^y^/a with 5n^y^ the 
usual Kronecker 5. Delta functions 5„.„^ in the periodic directions i > 1 are periodic 
of course. 



2.2.1 Quantization 

We now quantize the lattice Casimir system carefully and in detail, using the 
Schrodinger picture. On the lattice the classical field 0(x, t) becomes 

0(x,t) = 0(a;i,f_L,i) — > <Pn{t) = <j)n^,n^{t) ■ (2.27) 
subject to the boundary conditions 

<Po,nAt) = 07Vi,nx(t) =0, (2.28) 
0ni,nx(^) — 'Pni,n^+2Nj_ei(t) , i = 2, . . . , d . 

As before we ignore the external regions ni < and rii > Ni = l/a beyond the 
Dirichlet planes. Since = on these planes the interior lattice G on which 0(x) is 
a quantum field does not include these planes: 

G = {neZ'^\l<ni<{Ni-l);{-N^ + l)<ni<N^, i = 2,...d}. (2.29) 

From eq. ( |2.1| ) we obtain the classical Lagrangian 

^ = ^ E + E E ••• E ^ (0^ - 0n(m^ + n^i))0n - (V0)9 

^ [ni=0 ni=lj n2=-Af±+l na=-N^+l ^ 

= Y.^{^l-^n{m' + V{n{))c^^ + c^A/\4>)n) , (2.30) 
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where 

Afi-l A^x Nj_ 

E---=E E ■■■ E •■■ (2.31) 

neG ni=l n2=-Nj_+l na=-Nj_+l 

and A = An,n' is the lattice Laplace operator (see Appendix A). The specific 
ni-summation prescription given in ( p.30| ) avoids introducing an asymmetry in the 
Xi-direction by the asymmetry of the lattice derivative. The classical canonical 
momenta are 

^n{t) = -^ = a''Mt), ^neG (2.32) 

0(t)n{t) 



and the classical Hamiltonian is 

2 Ka?'^ 



H=Y.7Tn{t)Mt)-L=Y. ^ [^ + 0^(^2^V^(^^))0„_0„(A0)J .(2.33) 



nGG neG 

Canonical quantization in the Schrodinger picture assigns to the classical field (pn 
and momentum 7r„ time independent operators 

0„(t)-^$„, (2.34) 
a"- 

which satisfy the commutation relations 

= [n„,n„,] = o, (2.35) 

[^n,T^n'] = ^5n,n' ■ (2-36) 

Note that in the continuum limit a ^ 



^6n,n' ^ Six~x') (2.37) 



n„ ^ n(a;) . (2.38) 

The lattice Hamilton operator is 



n=J2^{ul + $„(m2 + \/(ni))$„ - <l'„(A$)„) . (2.39) 

ngG ^ 

We go now to the field representation in which the field operators are simply 
multiplicative: 

^nmM^t] = (2.40) 

Here \E'[{0m},^] is the lattice analogue of the wave functional \E'[0, t] of continuum 
theory, i.e. a function of the c-number variables 0m. The canonical-momentum 
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operators n„ become the differential operators — {i/a'^) d/d(j)n- In this representation 
the Hamilton operator assumes the form 

with (see eq. 

On,n' = ("^^ + ^(^l)) Sn,n' " A„,„/ 



d 



d 1 

V{ni))5n,n' -^ — {^^n+e,,n' -'^5n,n' + ^n-e,,n>) ■ (2.42) 



This matrix is symmetric in (n, n') and positive definite. We need its eigenvec- 
tors and eigenvalues to diagonahze the Hamilton operator. The eigenvectors {m^} 
fulfilling the eigenvalue equation 

Y: On,n'<, = elui . (2.43) 

n'sG 

can be written direct product: 

«n = <^n^■■< = <4l, (2.44) 

with 

k = {ki, /C2, • • • , kd) = {ki, k^_) . 

The indices k label the linearly independent vectors, of which there are exactly as 
many as there are points on the lattice - namely (A^^i — 1) (2A^^)'^~^ which is also 
the dimenson of the real symmetric matrix On,n'- One can therefore choose the 
individual fc's so they take the same values as the Ui and consequently form the 
same lattice: 

G = {keZ'^\l<ki<{Ni-l)]{-N^ + l)<ki<N^, i = 2,...d}. (2.45) 

Choosing the lattice as described above guarantees that the boundary conditions 
are automatically fulfilled. 

The separation formula ( ^.44| ) leads to (compare eqs. ( |2.10D - ( |2.12| )) 



5: [r(n05„,x - -(^„,+i,„;-2<5„,,„, +<5„,_iyJ)]t;^} = plvtl ,(2.46) 



n[=l 



J2 --^(5n,+l,n; - '25n,,n'^ + 5„,-l,n;) <5 = ^^^U^^^ , ^ = 2, ■ ■ ■ , (2.47) 



n'=-N^+l ^ 
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where , a;|. are related to the eigenvalues el of ( [^■43D by el = + pl + Yli=2 ■ 
While the can in general be obtained only numerically, the modes for the free 
directions can again be gotten analytically: 



2iVi 



4 
a 



(2.48) 



- sm 

2 ^2A^i 



In summary 



ki k\ 

V 



d-1 
2 



-fc I n I 



m 



Pk^+l^-^ sm 



i=2 



2Ni 



(2.49) 
(2.50) 



These formulae are the lattice analogue of eqs. ( p.l6| ), ( p.l7| ) in the continuum 
problem with one little difference: the exact lattice analog of the orthogonality- and 
completeness relation ( |2.18| ) would be Y^neG a-^u^ {i^nT = ^k,k' , Y^keG ""n i^n)* = 
8n,n'/of. In the limit a — > these expressions tend to (|2.18|) . But on the lattice it is 
more natural to work with the relations 



keG 



Sk,k' 
^n..n.' 



:2.5i) 

:2.52) 



It is therefore the combination u\l^~a^ which in the limit a — >■ tends to u^{x). 
To diagonalize the Hamilton operator we make the following expansion of 0„ in 
terms of the eigenvectors mJ^, 



keG V a 



:2.53) 



with the inversion 



E 

nGG 



Because the 0„ are real we have 



Eqs. (|2.53|) and ( p.54| ) imply for the derivatives 



d 



d 



E 

k&G 

E 



d 



d 



neG V a 



(2.54) 

(2.55) 

(2.56) 
(2.57) 
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In terms of the 0^ the Hamilton operator ( [^.411 ) becomes 

^ = ^ E (-T^f^ + 4 ^k4>l] , (2.58) 

a finite sum over uncoupled harmonic oscillator Hamiltonians. The remaining steps 
are familiar. We introduce annihilation/creation operators 

Ak = (e,,0,. + 4- I , (2.59) 



4 = ^ (^k^l - ,^-] , (2.60) 



with commutation relations [A^, AI,] = 6kk' , [Ak, A/.'] = [A^, AI,] = 0, in terms of 
which the Hamilton operator ( |2.58| ) becomes 

^ = ^ E (4^fc + AkAl) = Y.ek (AlAk + I) . (2.61) 

One recognizes that one is dealing with a finite set of uncoupled harmonic oscillators 
whose frequencies are the spectrum {e^}. A Fock space can now be constructed, 
whose ground state |0) is of course defined by |0) = 0\/k. The (finite) energy of 
this state is 

i^o = ^ E ^'^ • (2-62) 
The representation of (pn in terms of creation and annihilation operators is 

'^^ = E 4= f^^^ + ■ (2.63) 



2.2.2 Vacuum stress tensor 



a" V a 



To compute the vacuum expectation of the stress tensor operator ( p. 21 



we re- 



call that, in the Schrodinger picture and field representation, H(x,t) = (f){x,t) = 
—iTTr\ ■ Thus 

T'\x) = + (V0)^ + (m^ + ^(xi))0^) 

= \ [-^2 + + ^(^i) - + ^V(0V0) , (2.64) 



2 
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With the help of the correspondence (|2.38| ) and eqs. ( |A.32| ), (|A.33| ) one recasts 
( TO) , (EUD in the form 



1 92 



+ L — - 



n'eG 



i=i 



4a2 



(2.66) 



1 / 1 92 
2 



'n+Si ~ 



2^ c„,n'0n' - 2^ — 



n ' y^n—ei 



n'dG 



2 I ^2 



4a2 



2a2 



(2.67) 



The next step is to take the vacuum expectation values of the operators ( p.55|) , 
( |2.56|) . For this purpose it is useful to note that 



d 



d 



2£:fc 



(A 



A 



ki ,kj 



One then verifies 



(O|0,0*,|O) 



2ek 



d 



d 



|0) = --ekSk,k' 



(O|-^0fc,|O) = -(O|0fc-^|O) 



^k,k' 



(2.68) 
(2.69) 



(2.70) 



and 



r\2 

-— — 10) = {0\aUnj:On,n'M0) 



1 



U. 



k\2 



k€G 



(2.71) 



Using ( p.46| ), ( p.47| ) and the completeness relations for the eigenvectors one finds 

, (2.72) 



(oirio) = E^i4r+EE^(k..,r-2kr+k 
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pIi - ^(^i) , k 



? + y — = — 



u. 



-EE 

fceG j=2 



u. 



n+ei 



(2.73) 



and for i = 2, ■ ■ ■ , d 

mm = E 



UJ 



ki \„,k\2 



2a% 



+ E^(k...r-2kr+ 



fc |2 



[2.7A) 



Finally from eq. (|2.49|) 



\u 



k\2 



k^Y 



{2Ni_Y-^ 

so the preceding formulae can be simplified. Thus with L = a Nj_ 



(o|Trio) 

(0|T„i^|0) 



2ae, (2L)'^"i 



2(t;^J^ + «_i) 



8a%fc(2L) 



:2.75) 



^2.76) 



+ E 

k&G 



(oir|o) 



(^n^+l)^-2«)^ + (^;^,-l)^ 
8a%fc(2L)^-i 



(2.77) 



,5 2aefc (2L)'^-i 



- E 



«+i)'-2«)^ + (t;^,_i)^ 



8a3e,(2L) 



d-i 



i > 2 



(2.78) 



Eqs. (|2.76| ) - (|2.78|) are the lattice analogues of eqs. ( p.22|) - ( p.24| ) in the continuum 
problem. One can transform the latter directly to the lattice with the help of eqs. 
(|A.32|) , ( |A.33|) and the correspondence 



u' 



.,ki\2 
rii ) 



^2.79) 



Eqs. ( p.76| ) - (|2.78| ) are valid on the lattice n E G. However they do not hold for 
the Dirichlet boundary surfaces rii = 0, A^^i which have to be dealt with separately. 
On these surfaces one must use the unsymmetric derivative formula discussed in 
Appendix A which ensures a smooth continuation of the lattice mathematics onto 
the Dirichlet boundary surfaces. 
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The calculation of {T^'^) on the Dirichlet surfaces using the one-sided derivative 
is equivalent (see eqs. (|A.28|) , (|A.29|) ) to introducing in eqs. ( p.76|) - ( p.78|) the 
additional definitions 



[vt\Y ■■= , (2.80) 

-'A^i+l)^ (^A^i-l) 



:= (2.81) 



Thus on the Dirichlet boundary surfaces 

M\2 



(o|T°^,,jo) = (oinUlo) = E^^^^|^ 

We see that the unregularized (T^") and (T^^) are the same on the boundary surfaces. 
However, the renormalizations of (T^^) and (T^^) turn out to be different (sees. 3 
and 4). Thus the renormalized tensor components (T°°)ren and {T^^)ren on the 
boundary surfaces end up being different. 

Using the lattice definition ( |A.31|) of "spatial integration" along xi we find 



-■ / Ni Ni-1\ 

^ E + E E(0|dO) a'=j: ^e, = Eo . (2.83) 

^ \ni=l ni=0/ nx keG 

This relation must hold, of course, for (T°°) to be the unrenormalized lattice vacuum 
energy density. 

2.2.3 The case V{ni) = 

When the lattice potential V{ni) = between the Dirichlet walls rii = 0, A^i one 
can easily solve the lattice Schrodinger mode equation (|2.46|) . The eigenvectors and 
eigenvalues are 



<\ = xl M^hni) , (2.84) 



pI = lsin2(^), fci = 1,2,3,... ,iVi . (2.85) 



Thus 



k f2~ . Vrfci / 1 \ '2' i^fcxnj 



e 



2 2 4 . 2/^^1n 4 .2/ Tr/cj 

k ~ ~^ ' 
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The lattice vacuum energy is 



\ 



H — - sm 



7rA;i 



4 



sm 



^2.87) 



In computing the vacuum stress tensor we make use of the identities 

sin^(a (ra + 1)) - 2sin^(ara) + sin^(a - 1)) = 2 sin^(a) cos(2a ra) , (2.88) 



sin^(a) 



4 sin'' 



1 — sin' 



Thus 



sin2(!^(ni + 1)) - 2 sin2(^ni) + sin2(^(ni - 1)) 



2aVfc, cos 



2Tikini 



iVi 
2aVfci sin 



(2.89) 



2 2 • 2 / 



cos 



/ 2TTkini 



Then from eqs. ( p.76| ) - (|2.78| ) and with I = Ni/c 



(0|T°0|0) 



(o|r„"|o) 



/ {2Lr-^ 4^ 2 

\ - Pki 

k^ 21 {2LY-'e, 



sm 



2l{2LY-^ek 



cos 



27[kini 



iVi 



' Tcki 
2N, 



cos 



'27rA;ini 



^2.90) 



2/(2L)'^-i 



2 / Trki \ ( 2nkini \ 



{o\t:\o) 



V ^1 / 

1 V- 1 2 Pki + ^fc, (2-Kkin{ 



(2.91) 



Z (2L)'>'-i 2 



sm 



,2iVi. 



cos 



2'nkini 



iVi ; 

, z > 2. (2.92) 



Comparing eqs. (|2.90|) - ( p. 92 ) and the continuum formulae (|2.22|) - ( p. 24 ) with back- 
ground potential V{xi) = and v^^{xi) = ^j2/l sin^nkiXi/l), ki = 1,2,3,- ■■ one 
finds in the former the presence of a lattice artefact term 



nki 
^2Afi 



2'7rfcini ' 



(2.93) 



This term has no continuum counterpart and it diverges for a ^ 0. In Appendix 
B we investigate the ID version of A in considerable detail and conclude that it is 
essentially devoid of physical information. 
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3 One dimensional lattice without background po- 
tential 



In this and the following section we specialize to one-dimensional lattice systems 
with no additional free directions X2, ■ ■ ■ , Xd- We drop the subscript 1 and first 
consider the case of vanishing potential V{x). In sec. 4 we will introduce a specific 
nonzero potential along the lattice. In both sections we wish to deal with Casimir 
forces. This means we must also consider one or more regions external to the interval 
0<x</orO<n<A^ = //a. This is easily done. We can attach external regions 
which are exactly like the internal region < x < I (except for their length) with 
V{x) = and a Dirichlet boundary at either end. When we calculate the lattice 
vacuum energy E^aci}) in < x < / the same calculation gives us the lattice vacuum 
energy E^aci.L — 1) in the adjacent interval I < x < L. Likewise one calculation gives 
us the lattice T^^ in both intervals. Our principal goals here are to understand in 
detail how to renormalize E^ac and T^^ in the continuum limit a — ^ 0. Then we 
can calculate the Casimir force on the common Dirichlet boundary x = I m two 
ways: globally using E^ac and locally using T^^. Comparison with the continuum 
version of this system will be important. We shall find that in every respect the 
lattice theory closely approximates the continuum system as expected. It will be 
convenient to discuss the cases of zero and nonzero field mass separately. 



3.1 Renormalized vacuum energy 
3.1.1 Zero mass 

We find from eq. (p.87|) the lattice vacuum energy in the internal region Q < x <l: 

Ereg[h o) = - Pfc = - Sm ( — ) (3.1) 

^ fc=l k=l 

where = I /a. The label "res'" here (and subsequently) is used to denote lattice 
quantities. The finite sum (|3.1| ) can be evaluated using 

/, X fNx\ ({N-l)x\ fx\ , , 

2^ sm(fcx) = sm ) sm ( ) esc ( - ) (3.2) 



k=\ 

with the result 



4« 



Expanding in powers of a yields 
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Subtracting all terms in eq. ( ^.31 ) which diverge in the continuum limit a — we 
recover the continuum Casimir energy for a real ID massless scalar field confined 
between Dirichlet boundaries [3, Pi CT, 



lim 



21 



ira 



2 + 2a 



TC 



24/ 



(3.5) 



To directly derive E{1) from the lattice energy (|3.1| ) one can also do the following: 



lim 

AT— »oo 




.=1 v^sin(^; 



n \ s oo 
71 J 



Re(s) > 1 



(3.6) 



where ({s) is the Riemann C function |TT|. Because C(s) can be continued throughout 
the s-plane (and is a meromorphic function with its only pole at s = 1) one can use 
the preceding formula to evaluate the N oo limit of eq. (|3.1| ): 



1 Tik 
lim — - sin — - 

N-*ooaN 2N 



lim Ereg{l,a) 
N—^oo 



(3.7) 



Ivr 
72 



c(-i: 



71 

24/ 



To have a more complete physical system we add the external region I < x < L. 
Exactly the same calculations with / replaced by L — / give the lattice and continuum 
vacuum energies in this region; these are Ereg{L — I, a) and E{L — I) respectively. 
The total lattice and continuum vacuum energies in < x < L are then 



E{l,L,a) 



Ereg{l, o) + Ereg{L — /, o) 



(3.8) 



E(/,L) = E{l) + E{L-l) 



Similarly an external region —L' < x < could be attached to the other side of the 
interval [0,/]. Additional external regions beyond these could be attached. Global 
structure on x could also be imposed; e.g. periodicity on some large interval, or 
noncompactness for x — > ±oo. We shall do none of these things here. As stated, our 
sole purpose will be to probe the closeness of the lattice and continuum descriptions 
in < a: < L. 

Consider for example the global Casimir force on the Dirichlet boundary at x = /. 
The continuum definition of this force (in the direction of positive x) is 



d_ 



E{l,L) 



7T 



+ 



TT 



24/2 24(L-/)2 



(3.9) 
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On the lattice one would define 

Eregil + a,a) + EregiL — {I + a),a) — Ereg{l,a) — Ereg{L — l,a) 

i'Casyh 0.) = • 

a 

(3.10) 

In the limit a — ^ and using eq. (|3.3| ) one verifies that the same Casimir force (|3.9|) 
is recovered without the need to perform explicit subtractions. 



3.1.2 Nonzero mass 

For nonzero mass the lattice vacuum energy in the internal region < x < I 



1 ^"M 4 Tik 
EreJm,l,a) = i y Wm2 + i-sin2(— ^ 



1 irk 
= VM)^ + 4iV'sin2(^) (3-11) 

k=l ' 

evidently cannot be evaluated in closed form. It can, however, usefully be expanded 
in powers of ml: 

^ . , ^ N"^' ,7rk, 1 2,^1 1 

Ereg{m, I, a) = — }^ sm(— ) + -m l}^ 



I 8 iVsin(l^) 

/l/2\ r(3/2) 



z/!r(3/2-z/) 



The first term in eq. ( |3.12|) is just the m = energy ( |3.1|) which we already know 
how to renormalize. All of the terms in eq. (|3.12|) labeled by > 2 are convergent 
in the continuum limit. Each of them can, in this limit, be expressed in terms of 
the Riemann C, function using eq. (|3.6|) with s values s = 2v — 1 = 3, 5, ■ ■ ■ for 
which this latter series converges. Hence they require no renormalization. However 
the second term in eq. ( p.l2 ) does require renormalization. This is evident from 
eq. ( ^.61 ) where evaluation of the term in question corresponds to using ({s) at its 
pole s = 1. More directly, the second term in eq. (p.l2|) diverges logarithmically for 
— s> oo. Indeed for N ^ k 

, irk , nk , , 

iV=in(^)«- (3.13) 

and it is well known that 

lim fy i - IniV^ = 7 = 0.5772 . . . (3.14) 
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is a definition of Euler's constant 7. Analogously we define here 



N-l ^ 2 

Ett-T^ = - In iV + c(iV), (3.15) 
lim (c(A^)) = 0.52125... 

which specifies another sequence {c{N)} whose hmiting value for iV — >■ 00 can be 
obtained numerically with the result shown in eq. ( |3.15| ). Thus the second term on 
the right of eq. ( p.l2| ) behaves as 

-l.mHH^) + ^c{N) , (3.16) 

i.e. there is an additional logarithmic divergence which has to be subtracted. By 
subtracting from Ereg{m,l,a) the first term of ( ^.l(j| ) as it stands, we would incor- 
rectly discard the /-dependence of this term. We must therefore split this term into 
the finite term -^m^l ln(/i/) and a term —^rn?l ln(/ia) diverging for a ^ 0. Here 
is an arbitrary constant with dimension [mass] which has to be introduced to 
make the arguments of the logarithms dimensionless. We denote by E{m,l,a, ^) 
the expression obtained by subtracting this latter term, in addition to the divergent 
terms found in (|3.3|) , from Eregim, I, a), i.e. 

21 1 rn?l 

E{Tn,l,a,fx) = Eregirn,l,a) 9+7;- + -, — ln(/xa) . (3.17) 

Tea la 47r 

E(m,l,a, iJ,) converges in the limit a ^ 0. The renormalized continuum vacuum 
energy in the internal region < x < / can therefore be defined to be 

E(m,l, fi) := Urn E(m,l,a, fi) . (3.18) 

All the preceding global formulae can be transferred to the external region / < 
X < L by simply replacing I with L — I. The total vacuum energy for the system 
comprising the region < x < L with Dirichlet boundaries at x = 0, / and L is 
E{m, I, a, /i) + E{m, L — l,a, fi). The Casimir force on the boundary point at x = / 
is then 

Fcasi.'fn^l^L) = — lim - [E{m, I + a, a, fi) + E{m, L — I — a,a, ^) (3.19) 

— E{m,l,a, fi) — E{m, L — l,a, fi)] 

One can easily convince oneself that the divergences which have to be subtracted in 
order to make the renormalized continuum vacuum energy E{m,l,fi) finite, cancel 
in the difference on the right of eq. ( |3.19| ) in the limit a — > 0. I. e. in the lattice 



calculation of Fcas these terms do not need to be subtracted from Ereg and the mea- 
surable Casimir force, unlike the vacuum energy, does not depend on the arbitrary 
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parameter /x. 



Physically fixing /i 

The dependence of the m > renormalized continuum vacuum energy E{m,l,fi) 
on the arbitrary mass parameter /i seems unsatisfying. This does not, of course, 
mean that it is wrong. A vacuum energy in QFT evidently need not be absolutely 
calculable. (However shifts in vacuum energies - i.e. Casimir energies and their 
associated Casimir forces - should be absolutely calculable as is the case here.) We 
now propose a way to select a particular value of /i which, while not dictated by 
general theory, is physically motivated. 
To this end let us make the ansatz 

jj, = Km (3.20) 

where k is a dimensionless proportionality constant which will be determined by the 
requirement 

/ dE(m,L K,m)\ 



This means that the contribution to the Casimir force on the boundary a; = / is 

d_ 

dl 



determined by — ^ E(m, I, jj) alone if L is chosen to be oo. One can verify that for 
H = m 

lim E{m, I, n = m) = — — ln(x) - — , (3.22) 
i->oo 47r 4 

where x is a specific constant which can be determined numerically with the result 
Inx = —2.579436. Using this x we define a new renormalized continuum vacuum 
energy by 

E{m, I) = E{m, I, fj, = m) + — — In x 

( 21 1 m?l \ 

= Urn Ereg{m, I, a) i + 7^ + -r~ In(xma) , (3.23) 

y 7ra^ za Air J 

which has the large-Z behavior 

TTl 

Mm E{m,l) = -— . (3.24) 

Hence, if we identify the parameter k introduced above with the numerically deter- 
mined constant x? this vacuum energy E{m, I) achieves the goal set in eq. ( p.21 
Let us give a bit more detail. From eqs. (|3.5|), ( p. Ill ) - ( |3.15| ) 



n 1 

E{m,l) = -— + —m^l\n{xml) 
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Figure 1: E{m, I) from (|3.23|) with x = 0.0758166. E{m, I) is summed to 
k = 1000. 

In this and the following figures an arbitrary length unit is used. 



+ 



— lim 

21 N^oo 




, irk , 



-2A^sini 



vr/c 



[mlf 



'2N' ANsini 



2N' 



IT 

24/ 



+ j-mH Hxml) + 1 f; (^^{miy + (nky - Trk - M!^ (3.25) 

where ln(x) differs from ln(x) by 7rc(oo)/2, c(oo) being the constant introduced in 
(EH): 



ln(x) =ln(x) + - 0.52125... 



-1.76066... . 



(3.26) 



In fig. 1 we plot i?(m, Z) versus / for three different masses. An arbitrary length 
unit is involved in Fig. 1. The relevant scaling property of E{m,l) is obvious in 
( |3.25| ). In natural units {h = c = 1) the only dimensional quantities available are I 



and m, with ml dimensionless. Thus I sets the scale in ( |3.25| ). E(m,l) is negative, 
increasingly so for increasing m and fixed /. For fixed m and increasing / we see 
that E{m, I) increases, corresponding to an attractive Casimir force Fcas between 
the boundaries at x = and x = I, where due to our fixing of /x no account has to 
be taken of the external regions x < 0, a; > / provided they are chosen to extend to 
infinity. In this case the Casimir force on the boundaries at a; = 0, / becomes 



N~l 



Cas 



'm,l) 



lim 



sm 



2( -Kk^ 
-2N ' 



k=l 



2/ Trfc > 
2N' 



1 

2d 



2 2 \ 

-Hxma)--j 



(3.27) 
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Comparison with analytic results 

Many authors have calculated the vacuum energy shift of the massless and massive 
ID scalar field confined between Dirichlet points. An early and useful analysis in 
arbitrary dimension was given by Ambj0rn and Wolfram 0. We simply quote their 
result for d=l: 

EAw{m = OJ) = (3.28) 
EAw{m,l) = -TT 1^ 



lim EAwi^y = , 



l—fOO 



where Ki{z) is a modified Bessel function. Note that no arbitrary mass param- 
eter fi appears in Eaw{'^J) and that limi^oo E aw J) i-e. the situation is 
essentially the same as above after fixing ^ by the requirement ^ i?(m, /) = 0, 
the only difference being that our E{m,l) tends to — m/4 for / — oo whereas 
lim;^oo EAw{fn, = 0- Let us form the ratio 

EAwim,l) 12ml ^ K,{2mln) 

eAw{ml) = — = — — 2^ 3.29 

EAw{m = 0,l) TT^ ^1 n 

and the corresponding ratio using our result E{m, I): 

_ EimJ) + ^ Qml Qjmlf , 
'^"^^^ = Eim = OJ) = 

- v£fvWTW--A:-M!) . (3.30) 



^ fc=i 



Numerically very close agreement of these ratios up to (9(10^^) - (9(10~^) is found. 
Their difference increases for large ml as (ml)"^; this arises from our imperfect knowl- 
edge of the number In^. 

3.2 Renormalized vacuum stress tensor 
3.2.1 (T") 

From eq. (|2.91| ) we find for the internal region 

( Til ^ ( i^ri,a) = -Y: I ^^^^ , , + A{m, /, n, a) (3.31) 
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where the lattice artefact term 



A{mJ,n,a) = V] , ^ ^ ' (3.32) 

ti J(m/)2 + 4iV2sm2(^) 



has already been mentioned. This term has no counterpart in the continuum theory 
and it is unphysical. As discussed in Appendix B it makes no contribution to the 
Casimir force. The artefact A diverges as a — >• and must be subtracted - but not 
yet for there is some interplay (or more precisely a cancellation) between the leading 
divergence in A and a nonleading divergence in the first sum in eq. ( p. 311 ). 



For zero mass m = one can express eq. ( |3.31| ) using eqs. ( |3.3|) , ( p.l8|) in the form 
( Til ^^^^(^^ = ^ _ _^ + a) + 0{a') , (3.33) 

where the lattice function f{l,n,a) is defined by eq. ( [B.19| ). Because / ^ as 
a ^ in the internal region < x < I the renormalized continuum tensor for m = 
in this region is 

as it should be (recall eq. ( p.9| )). For any m > the continuum (TH) is constant 
between parallel Dirichlet planes in arbitrary dimension (see e.g. ref. fl^). This 
constant value of (TH) equals the Casimir force/area on the enclosing Dirichlet 
boundaries when no force acts from outside on these boundaries (because e.g. the 
external region is infinite). 

At the lattice level note that the physical first term in eq. ( |3.31| ) is indeed indepen- 
dent of lattice position n. Thus the lattice formulation has already scored a success 
in reproducing this known continuum feature: 

2 ^"1 sin^f^l 

(r" W = (r" = E 7=— (3.35) 

" fc=i ^(m/)^ + 4A/^ sm^(l^) 

On the other hand the lattice artefact second term in eq. (p.31|) does depend on n 
- another indication of its unphysical nature. 

It is not difficult to see that {T^^){m, I, a) coincides with the nonrenormalized Casi- 
mir force obtained from the nonrenormalized global vacuum energy. Thus on the 
lattice one obtains the same Casimir force using global and local methods, just as 
one does in the continuum formulation. Because we know the a — >■ divergent 
structure of eq. ( p.35| ) from the global problem we can write down immediately the 



renormalized result 

(Tii)(m,/,/i) = ^ J: ^===::^^ML=== (3.36) 




+ — - — ln(/ia) 



TTa^ 2al An 
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Of course the arbitrary parameter /i is again present. 

Physically fixing /i 

The condition chosen here to fix fi is 

{r")KO = -^^ (3.37) 

where E{m,l) is the fixed- vacuum energy ( |3.25|) . Equivalently 

lim (T")(m,/) = . (3.38) 

It is important to emphasize that in general 

{T''){m,l,fi) ^ , (3.39) 

assuming one uses the same /i on both sides. Directly from eq. (|3.26| ) we obtain the 
fixed-/i vacuum tensor 

(r")(™.o = -%i> (3.40) 

— lim I — \ — ^ ^ — -|- — — \]i(t7ti(i) 

^1 ^{miy + Am sm\^) ^ra^ 2al An 

where 

ln(r) = ln(x) + 1 = -1.579438 .... (3.41) 

Here fixing /i has led to fi = rm. In the global vacuum energy it was n = 

With the V{n) ^ system of sec. 4 in mind let us say something more about 
the renormalization of eq. ( p.31| ). Subtract the lattice artefact term A (using eq. 
( |B.26| ) for it) in addition to the other subtractions in eq. ( |3.36| ): 

/ 2 m} 
l,T^^){m,l,^) = hm (T") (m,/,n,a) ^ - ln(/ia) 



Tn?' I \ 
-/(/,n,a)-— /i(-,n) , (3.42) 
8a/ 



or with n = rm 



( 2 

'T^^){m,l) = lim (T^i) (m,/,n,a) - — In(rma) 

a— >o y y vra^ Att 

-f{l,n,a)-—hi-,n)) , (3.43) 
8a/ 
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where /, h are defined in Appendix B. As indicated both of these functions depend 
on position x = na. However, in the continuum hmit they both vanish everywhere 
except on the boundaries n = 0, N where they take values independent of /: / — 
const /a? while /i — > 2/7r. Subtracting this position dependence in eq. ( |3.43| ) yields 
a (T^^) independent of lattice position. To show this we rewrite eq. (|3.43|) in the 
form 



N-l 



{T'^)(m,l) = limV^ ^=^^M^_ - 11 sin(^) + ^--^ 



^ ^ ^ i / _l sin f ) cosf ) 

/y(m/)2 + 4Ar2sin2(|±) I' '2iV^ / 

,TTk, ,2TTkxA TT m? 



sin(— ) cos(^) - -— - — ln(fm/) (3.44) 




24/2 2/2 \ ^ l(^]\2 j_ 2'nk 



with 



ln(f) = ln(r) + ^0.52125... (3.45) 
= ln(x) + l = -0.76066... . 

The two terms following the first term in curly brackets are the first two terms of 
its expansion in m? multiplied by (-1). This removes all a — > divergences and the 
resulting limiting value (second equality) follows straightforwardly. Expansion of 
the latter in powers of ml yields 

(T")(m,/) = -_ln(fm/) + — ^rfJ— j Ci^v-l) , with (3.46) 
d. ^ (_,).r(l/2 + .) 



V V\ \/TT 



In fig. 2 we plot the ratio {T^^){m,l)/ {T^^){m = 0,1) versus ml. The exponential 
fall-off of this function with increasing ml is evident. 



3.2.2 (TOO) 

From eq. (|2.9C ) we have 



(TOO),,^(m,/,n,a) 



Ereg{m, I, a) m2 ^I-^^ 

/ 2~ 

^ ^ k=i 

+A{m, I, n, a) . 



cos(^) 



^{miy + AN^ sin^dl) 



(3.47) 
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Figure 2: The ratio {T^^){m,l)/ {T^^){0,1) with Inr = —0.76066 versus ml. 
Eq. (|3^) is summed to k = 1000. 



The property 



n=l 



for A; 7^ 



(3.48) 



then yields immediately (compare eq. ( |2.83|) ) 



N 



E ( )reg("^' ^' ^' = ^reg(m, /, o) . 



(3.49) 



n=l 



Let us now renormalize eq. (|3.47|) . 

The second term in eq. (|3.471 ) remains finite in the limit N —>■ oo away from the 
boundaries 



'N-1 

lim 



COS I 



E 



cos 



• 2-Kkx ' 



fc=i ^(m/)2 + 4iV2sin2(i|); ^tl y'(m/)2 + {j,kf ' 



(3.50) 



This series converges for < x < /. However its behavior near the boundaries 
is Inx and ln(Z — x) for a; ^ and x — > I respectively (where the series diverges 
logarithmically). The renormalized energy density 



^ 2 y^(m/)2 + (7rA;)2 



(3.51) 



Tra^ 47r 



a a 
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Figure 3: The function eq. ( p.54[ ) summed to k = 1000 over the interval 

0<x/l < 1. 



therefore contains the expected boundary divergences. Fixing /i by = xm we 
obtain 

,ooe.,U)^^-^£-J2lttL, (3.52) 

= f ( ),,^(m, lA,a)-—^ + ^ In(xma) - /(/, ^, a) - ^M^, . 
a^o \^ '^y a Tia'^ An a S a a J 

Summing eq. (|3.52|) over the lattice < x < / we find with the help of eqs. ( p.ll|) , 



H « ( )reni^^ i,n,a) 2 ^ IT Mx^tt) - f{l, n, a) - —h{-, n) 



Ereg{m,l,a) - — + — + -^\n{xma) , (3.53) 



n=l 

21 1 mH 

which is the renormalized global vacuum energy ( p.23|) as it should be. This entitles 
us to assign parameter x the same value it has in eq. (|3.23|) . Here we do not have 
a condition to impose which fixes x do we need one. The connection (|3.53|) is 
quite sufficient. 

The energy density (T°°) is in contrast with the constant (T^^) strongly position 
dependent, especially important being logarithmic divergences at the boundaries. 
To display this we plot in fig. 3 the quantity 

00. / 2; _ (ro°)(m,/,x) E{m,l) 
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12 (ml)^ ^ cos(gf^) 

y^(m/)2 + (7rA;)2 

versus x/l for three values of ml. The boundary divergences in (T°°) are inseparable 
from the boundaries themselves [Q. It is gratifying that the lattice formalism very 
closely reproduces 0, ^ this crucially-important and prominent feature of contin- 
uum QFT. 



4 One dimensional lattice with Bessel potential 
background 

Continuing on the ID lattice with Dirichlet point boundaries at xi = 0, X2 = I 
and xs = L (in this section it will be convenient to make this change of notation) 
we now introduce a background potential V{x) > 0. The Dirichlet boundaries at 
X = Xi, X2, X3 are embedded in V{x). For convenience we shall in this section 
combine V{x) and the mass term into a single potential 

U{x) = m^ + V{x) . (4.1) 

In subsection 4.1 we compute the lattice vacuum stress tensor and vacuum energy 
for the region Xi < x < X2 between adjacent Dirichlet boundaries embedded in an 
arbitrary background potential V{x) which vanishes as x — > 00. Then in subsection 
4.2 we choose V{x) to be the potential 

V{x) = % , a > 1/2 (4.2) 

{x - Xo)^ 

and begin the principal numerical work of this paper. We compute numerically 
(T^'^) in xi < X < X3 and the global vacuum energy in this same region. Then we 
compute the Casimir force on the middle Dirichlet boundary at x = X2 = I first as 
the jump in (T^^) at this point, and second using the global method. Very close 
agreement between these two forces (which should be the same, of course) is found. 
This is a crucial test of our lattice method, and it succeeds. A variety of other checks 
and tests will also be displayed in subsection 4.2. 

We remind the reader that the potential ( |4.2|) has a specific physical interpreta- 
tion. It defines what the authors call a "Bessel boundary" ^ - i.e. semitrans- 
parent surface texture represented by V{x) attached to a core Dirichlet boundary 
at X = Xq. This entire boundary was absent in the analysis of sec. 3. Parameter 
a > 1 /2 controls the amount of semitransparent material present and the potential 
V{x) extends from xq into x > xq as far as we choose to let it. [We can also let V{x) 
extend arbitrarily far into x < xq - and in fact should do so if we wanted to calculate 
the Casimir force on the Bessel wall - but so far as this section is concerned, we 
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need not bother about the region x < xq-] Without the quantum field the classical 
background material represented by V{x) would not interact directly with the clas- 
sical boundaries at x = xi, X2, x^ nor would these objects interact with one another. 
However, when the quantum field is present it is distorted by V{x) and by each 
of the hard boundaries at x = Xi, X2, x^. Back forces are then exerted by on the 
extended structure V{x) and on each of the hard boundaries. Forces involving all 
parts of the spatial background (including the material represented by V{x)) come 
into existence. Thus we have a highly nontrivial Casimir system to study. 



4.1 Renormalized vacuum stress tensor for semihard V{x) 
4.1.1 (T") 



For the region Xi < x < X2 we have from eqs. (|2.77| ), the lattice formula for arbitrary 
background potential V{n) 

{T'')regi^i,X2,x^ + na,a) = \\ (4-3) 

fe=l \ ^"-^fc 

fci)'-2(t;^)^ + (^;ti)^\ 



+ 



Here X2 — X1 = Na = I. As we show in Appendix C this expression can be renormal- 
ized for arbitrary V{n). The renormalized lattice (T^^) is (prior to letting a 0) 

2 

{T^'^){xi,X2,Xi + na,a,fi) = ( ) (xi, X2, Xi + na, a) - f{l,n,a) 

U (xi + na) , , , U(xi+ na) . ,1 . , , 
^ ' ln(/ia) ^ /;,(-, n) (4.4) 



47r ' 8 a 



where h and / are defined by eqs. ( |B.12| ), ( |B.21D respectively. As in the corre- 



sponding lattice formula (|3.36|) for Vin) = the renormalization process introduces 
dependence on an arbitrary renormalization mass parameter fi (as always because 
of the need to subtract a lattice sum diverging logarithmically as the lattice con- 
stant a — ^ 0). Such an arbitrary renormalization parameter introduces ambiguity, 
of course. Fortunately, because of the way (T^^) depends on fi in eq. ( [4.4| ) this pa- 
rameter does not contribute to Casimir forces on e.g. the boundaries at x = xi and 
X2- The latter are determined by the jump in (T^^) across these points. If the back- 
ground potential ^(x) is continuous across a boundary, then the term U{x) In /ia 
takes on the same value on both sides of this boundary: there is no jump in ( T^^ ) 
and therefore no Casimir force. Only if V{x) makes a jump at a boundary will 
the Casimir force on this boundary depend on fi. Because Casimir forces are the 
principal experimental signature of the distortion of the quantum field by its inter- 
action with all background structure the importance of Fcas not depending on /x 
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for continuous background potential V{x) is clear. This observation also highlights 
the importance of not merely considering the internal region in studying a Casimir 
problem. To have a well-posed problem external regions have to be considered as 
well. 

Quite another question to be answered is this: Is the renormalized vacuum stress 
tensor itself well-defined? {T^'^) is after all our basic local mathematical instrument 
for representing/investigating the physical properties of the system. If n in eq. (|4.4|) 
is truly arbitrary then {T'^'^) is not unique. What are the consequences of this? Is the 
nonuniqueness associated with /i relatively trivial or is it symptomatic of something 
deeper? Even if does not influence Casimir forces it does influence the energy 
density (T°°) and the total vacuum energy. 

We see no compelling reason to view n as being, in some sense, a universal con- 
stant. If it were, then why would it be absent for U{x) = 0? Thus in sec. 3 where 
V{x) = we have fixed this ambiguity by an ad hoc but reasonable condition 

lim (T^^)(m,/) = . (4.5) 

Physically this means that there is no contribution to the Casimir force on the hard 
boundary at x = Z coming from the infinitely extended region x > I. 
Here with x-dependent potential V{x) > we choose what we believe is the natural 
extension of the scheme in sec. 3. We set 

= r^u{x) = [m^ + V{x)] (4.6) 



which of course reduces to the ansatz fx = rm for V{x) = 0. Then eq. ( ^4.4|) becomes 

2 

(T^^)(xi,X2,xi +na,a) = ( T^^ ) (xi, X2, xi + na, a) - f{l,n,a) 

U{xi +na) 2tt/ . \ 2\ U{xi+ no) , J . , , 

^ ' \w{r^\J [xx ^ na)a^) ^-^ -h{-,n). (4.7) 



87r 8 a 

Note that r here still does not influence any hard-boundary Casimir force as long as 
V{x) is continuous. Because r is a constant we can fix it by going to spatial infinity 
where (we assume) V{x) ^ so that U{x) m^: 

lim lim (T")(xi, xs, x, u) = . (4.8) 

Here Xi < x < X2 prior to moving X2 = / out to infinity and then following with x. 
This condition is essentially the same as that posed in the V{x) = case in sec. 3 
and, just as in eq. ( p.41| ), we obtain the fixed r 

Inr = -0.76066 • • • - | 0.52125 . . . = -1.579438 .... (4.9) 

Eq. ( [4.7| ) generalizes eq. ( |3.43| ) in the sense U{x). Now in the continuum 

limit with x = xi + na held constant, / = X2 — xi and x = x — xi 

(Ti^)(xi,X2,x) = lim((Tii) (xi,X2,x,a)-^-/(/,-,a) (4.10) 
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ln(r"'f/(x)a^) —h{-,-)] . 



Stt S a a ^ 

4.1.2 (TOO) 

The lattice formula for (T°°) is 

The renormalized lattice vacuum energy generalizing eq. (|3.53|) is 

2 

(T°°)(xi,X2,Xi + na,a, /i) = ( T°° ) (xi, X2, Xi + na, a) --f{l,n,a) 

^ln ;,V '-h{-,n). 4.12 

Stt 8 a 

With = ?7(x) 

2 

(T°°)(xi,X2,Xi + na,a) = ( T°° ) (xi, Xs, Xi + na, a) - f{l,n,a) 

U(xi + na) ^ , , r,. U(xi+na),.l 

Stt 8 a 

where 

In X = Inr - 1 = -2.579438. .. . (4.14) 

In the continuum limit 

(T00)(xi,X2,x) = lim('(T00),,^(xi,X2,x,a)-^-/(/,-,a) (4.15) 

+ ^ Hx'U X a^ _ ^ - . 
87r Q a a J 

4.1.3 Renormalized vacuum energy 

The lattice vacuum energy in the interval xi < x < X2 

^ N-l 

EregiXi,X2,a) = - J2 (4.16) 
^ fc=l 

is renormalized in Appendix C with the result 

£'(xi,X2,a,/i) = i?reg(xi,X2, a) - ^^^^ + tt" + / U{x)dx ln(/ia) .(4.17) 

vra"' 2a 47r Jxi 
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Assuming E = J {T^^)dx we find from eq. ( [4. 131 ) tlie same renormalized lattice 
vacuum energy 



E{xi,X2,a) = Ereg{Xl,X2,a) - 
1 f^^ 



2{X2-Xi) 1 



7ra^ 



+ 7- r U{x)\n{x^U{x)a^)dx . (4.18) 

OTT Jxi 



Note tliat in eqs. ( |4.17D , ( [4. 18]) one can use in place of the integral the lattice sum 

1 /N-l N \ 



IGtt 



^ + ^ a f/(a;i + na) ln.{x^U {x I + n(i)a ) 



(4.19) 



(4.20) 



^?i=0 n=l/ 

The difference vanishes when a 0. In the continuum limit we define 

E(xi, X2) := lim E(xi,X2, a) . 

4.2 Numerical evaluation using the Bessel potential 

Now we proceed to a numerical investigation of the renormalized (T^'^) and E{xi, X2) 
defined above for the background Bessel potential (|4.2|) . The core of the Bessel 
boundary is at x = xq and the surface texture represented by V{x) extends to the 
right from x = xq- We generally choose xq = —0.01 and position the first Dirichlet 
boundary at Xi = 0. The other two Dirichlet boundaries are further to the right 
ai X = X2, X3 > 0. Before presenting our numerical results it is perhaps worth 
saying again how they are obtained. First the lattice mode equation ( |2.46|) has to 
be solved numerically for the lattice modes and their energy eigenvalues Van- 
ishing boundary conditions at each of the hard boundaries are additionally imposed 
on the v^. A given calculation is performed for specific values of lattice constant 
a, of parameter a in eq. ( [4.2|) , of mass m and of positions Xi,X2,X3. When any of 
these six values is changed the numerical calculation has to be redone. Given the 
numerical ingredients {f^} and {e^} one can insert these ingredients into eqs. ( [4.4|) , 
([4.13|) to obtain the lattice (T^'^) numerically. 



4.2.1 Renormalized (T°°) 

Fig. 4 shows the function (T°°) obtained from eq. ( |4.13| ) in xi < x < X2 with 
Xq = —0.01, Xi = 0, X2 = 2 for m = and three different values of the parameter 
a in eq. ([4.2| ): these are a = 1/2 (or V{x) = everywhere between the Dirichlet 
boundaries), a = 1 and a = 5. Note again the arbitrary length unit along the 
horizontal axis. For a = 1/2 we recover the well-known (m = 0)-result (see e.g. 
ref.jig) (r°°) = 7rV(24/2). For a > 1/2 we see that (r°°) is (i) asymmetric 
about the midpoint x = 1 because of the asymmetry of the potential ( [4.2| ) and 
(ii) divergent as the boundaries are approached. To investigate this divergence 
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Figure 4: (T°°) from eq. ( [4.13|) with m = 0, Xq = —0.01 for Dirichlet bound- 
aries at xi = and X2 = 2. The lattice constant is a = 1/400. 
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Figure 5: (T^^) from eq. ( |4.13| ) with a = 1 and other parameters as in fig. 4, 
evaluated on the right-hand boundary x = X2 = 2. 
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Figure 6: (T^^) from eq. (|4.7| ) with xq = —0.01, a = 1 for Dirichlet bound- 
aries at Xi = and X2 = 6. The lattice constant is a = 1/300. 



differently we display in fig. 5 (T^^) for n = or x = X2 = 2 on the right hand 
boundary and a = 1 as a function of (—In a). Clearly a logarithmic divergence is 
involved. 



4.2.2 Renormalized (T") 



Here we use eq. ( [4.7| ) with the potential ( |4.2| ) to compute numerically the renor- 
malized (T^^). First of all we wish to show what this function looks like when the 
Dirichlet boundaries are absent and the system consists of the scalar field coexisting 
with the Bessel boundary. Of course this system can only approximately be realized 
on a finite lattice. Again we position a Dirichlet boundary at x = xi = very close 
to the core xq = —0.01 of the Bessel boundary. This is done for computational 
reasons. We expect (and we may also conclude from our numerical results) that this 
hard boundary essentially disappears into the (Dirichlet like) Bessel boundary and 
has no significant effect on (T^^) away from x > 0. Also for computational reasons 
we introduce a second Dirichlet boundary at X2 = 6 "far" from x = 0. We obtain 
the numerical results shown in fig. 6 for three different mass values m = 0, 2, 5 (and 
a = 1, a = 1/300). These need to be discussed in some detail. 



IS 



For "large" distance x from the core of the Bessel boundary the potential ( 4.2 
essentially constant (and small). Thus one expects (T^^) for large x to approach a 
constant value equal to (T^^) in free space for the same field mass m (see sec. 4.2.1). 
This behavior is evident in fig. 6 for masses m = 2, 5. For smaller mass m = 2 the 
asymptotic value for large x is close to zero. For larger mass m = 5 the asymptotic 
value of (T^^) increases noticeably. We remind the reader of the free-space vacuum 
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stress tensor in one spatial dimension (see e.g. ref. ||1U||) 



T;;,, = r^^^A, A = -— In- . (4.21) 
■' 47r /i 

This formula indicates that, if one may assume m > /i, Tjl^^ is positive and increas- 
ing with increasing m - just what is observed in fig. 6. But remember that we fixed 
r in eq. ( |4.7D by requiring (T^^) to vanish at large distance. (The same requirement 
applied to Tj^'ee would lead to fi = m and this tensor vanishing everywhere.) This 
is observed in fig. 6 for m = 2 but not for m = 5. The mathematical reason for the 
m = 5 curve not vanishing is our imperfect knowledge of the value of parameter r. 
The numerical consequences of this grow with m, one consequence being that (T^^) 
does not quite vanish for large x although it should. Note also that with increasing 
m the asymptotic value of (T^^) is reached in a shorter x interval. This is entirely 
consistent with continuum results where boundary effects are known to fall off ex- 
ponentially in mx. 

The (m = 0)-curve in fig. 6 does not approach the value (T^^) = asymptotically, 
nor should it. For m = there is no parameter to fix and no asymptotic condition 
on (T^^). For m = boundary effects in (T^^) no longer diminish exponentially 
away from the a; = boundary, but rather diminish as the inverse power Con- 
sequently the boundary at X3 = 6 is no longer "far away" and in fact is strongly 
influencing (T^^) in fig. 6. To see (T^^) asymptotically approaching zero one must 
choose a substantially larger value of X2- The m = continuum vacuum stress 
tensor between Dirichlet boundaries at x = 0, L is 

= -5,. ^ . (4.22) 

In fig. 6 the (m = 0)-curve for (T^^) goes negative as it should. 
Now let us modify the system by introducing a third boundary ai x = X2 = 1, 
leaving the other two boundaries in place at xi = and X3 = 6 (the previous X2 is 
now renamed X3). Figs. 7, 8 display (T^^) throughout the interval < x < 6. One 
sees in fig. 7 that Fcas on the intermediate boundary X2 = I defined by 

is directed to the left. This is predictable: the Bessel boundary's core Xq = —0.01 
and the Dirichlet point Xi = are much closer to X2 = 1 than is the other Dirichlet 
boundary X3 = 6. Dirichlet objects attract in Casimir theory and the Casimir force 
between closer Dirichlet objects dominates. Also predictable is the decrease shown 
in fig. 7 of Fcas with increasing mass m. For V{x) = between Dirichlet boundaries 
it is known that Casimir forces weaken rapidly (exponentially) with increasing mass. 
There is no reason why any background potential V{x) of the type considered here 
would qualitatively alter this general rule. 

Fig. 8 shows (T^^) for fixed m = and three values of the strength parameter 
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Figure 7: (T^^) from eq. with xq = -0.01, a = 1, a = 1/300 for Dirich- 

let boundaries at Xi = 0, 0:2 = 1 and X3 = 6. 
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Figure 8: (T^^) from eq. ( [4. 71) with xq = —0.01, m = 0, a 
Dirichlet boundaries at xi = 0, X2 = 1 and 0:3 = 6. 
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a = 0.5, 0.6 and 1.0 in potential V(x). Clearly Fcas increases with increasing a. This 
too is predictable. When a = 1/2 and V{x) = we know that (T^^) = -Tr/{2AL^) 
is constant (with L = 1 to the left of X2 = 1 and L = 5 to the right of X2 = 1). 
Increasing a above a = 1/2 causes V{x) =0 to jump abruptly and nonuniformly 
to a positive value everywhere, and obviously to a much larger average value in 
< X < 1 than in 1 < x < 6. The potential V{x) represents "Dirichlet material" 
added throughout the interval < x < 6. This nonuniformly added material pulls 
the X2 = 1 boundary more strongly to the left than to the right. 



4.2.3 Re normalized vacuum energy 

These numerical calculations of the vacuum energy are based on eq. ( [4.18|) , where 
the integral can be done analytically for the potential ( [4.2p . With I = X2 — Xi and 
6 = — 1/4 we find 



X2 



U{x) ln{x^U{x)a'^)dx 



In 



P 



X2 



U{x) dx 



(4.23) 



+ - / U{xi + xl) f ln([/(xi + xl) f) dx 
I Jo 



and 



U{xi,X2) 

g{xi,X2) 



1 f^^ 
1 



y (x - Xo)2 

,.2 1 



dx = m + 



a2 - 1 

" 4 



[Xi - Xo){x2 - Xo) 
,2 1 \ 



(4.24) 



dx 



{Xi - Xo)ix2 - Xo) 



+Aml\/b 



arctan 



'm {x2 - Xo)' 



— arctan 



' m (xi — Xo) 



I {m?{x2 - XqY ~^\^( ^'^{^^{^2 - xqY + hY 



X2 - Xo 



1X2 - Xo 



I {m?ixi — Xo)^ — ^) , / /^(m^(xi — Xq)^ + h) ~ 
In - 



Xl - Xo 



(xi - Xo)^ 



The renormalized vacuum energy in region Xi < x < X2 is 



N-l 



£^(xi,X2,a) = -^Sk 



k=l 



21 1 
7ra^ 2a 



lU{xi,X2) ^^^Xa^ ^ g{xi,X2) ^^25) 



Air 



I 



8nl 



Fig. 9 displays this energy as a function of boundary position X2 for fixed boundary 
position Xl = 0. Nothing beyond the boundary at X2 is considered here. Fig. 9 
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Figure 9: E{xi, X2, a) from eq. ( [4.25|) with Xq = —0.01, m 
for Dirichlet boundaries at xi = and variable X2- 



2, a = 1/500 



qualitatively resembles fig. 1 giving the comparable vacuum energy for V{x) = 0. 
As X2 approaches X2 = here (the boundaries Xi = and X2 move close together) 
the function E{0,X2,a) approaches the function — 7r/24/ with I = X2 — xi. This is 
the / — > limiting form of the vacuum energy for V{x) = and m > as it should 
be. 



4.2.4 Casimir forces 

Now we proceed to the numerical evaluation of Casimir forces for various configura- 
tions. It will be extremely important to verify that we obtain locally from (T^^) and 
globally from the vacuum energy the same Casimir forces. If this is what is found - 
and it is - then because the local and global mathematics differ quite substantially 
we can be confident both are working as they should. Indeed these calculations 
provide the crucial test of our lattice methods with nonzero background potential 
V{x). 

Global calculations of the vacuum energy will be done using the renormalized for- 



mula ( [4.25| ). Positioning three Dirichlet boundaries as before at Xi < X2 < x^ we 
obtain the Casimir force on boundary X2 (globally) as 

E{xi, X2, a) — E{xi, X2 + a,a) + E{x2, X3, a) — E{x2 — a, X3, a) \ 



Fcas = lim 



Alternatively Fcas can be calculated (locally) as 



(4.26) 



Fcas = \im{{T^^)^^g{xi,X2,X2,a) - {T^^)^.^g{x2,X3,X2,a)) . (4.27) 
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Figure 10: Casimir force on the middle boundary X2 = 1 obtained globally and 
locally from eqs. j ^M ), (1^ ) with xq = -0.01, m = 0,a = 1/300 



and xi = 0, X3 = 6. 

Fig. 10 shows these two results for Fcas for fixed xi,X2,X3 as a function of Bessel 
parameter a. The agreement is very good. Let us examine this in greater detail. 
For small a we find very close agreement between the local and global calculations 
even for relatively large lattice constant a. For example, in fig. 10 at a = 1 the 
global and local values are -0.246007 and -0.246225 respectively. The limiting value 
of Fcas for a -> 1/2 is -7r/24 + 7r/600. 

As a increases the agreement between global and local becomes less good. To im- 
prove the agreement one must go to smaller lattice constant a. This is demonstrated 
in fig. 11 where, for a = 10, we plot the globally and locally obtained results for Fcas 
versus a. Linear extrapolation of the two straight lines to a = yields the limiting 
values -2.45142 and -2.45105 respectively. This is extremely good agreement. 

Note that in fig. 10 the (m = 0)-Casimir force is seen to depend linearly on a. 
For m > this is not the case in general; however linear behavior does set in for 
sufficiently large a. 

To display the agreement between the global and local methods in even more detail 
let us decompose Fcas on the X2-boundary into left (L) and right (R) components 

Fcas = Fl + Fr . (4.28) 

Fl{Fr) is the force exerted on the boundary X2 from left (right) 

Fl ^ -^^^^^ = {T''){x„X2,x = X2) , (4.29) 
0x2 

_ dE{x2,X3) /rj.u\f \ 

Fr = ^ = -{T''}{X2,X3,X = X2) . 

0x2 
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Figure 11: The global and local Casimir forces in fig. 10 for a = 10 plotted 
versus lattice constant a. 
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Figure 12: The left-hand Casimir forces ( 4.29| ), ( [4.30| ) on the middle wall, the 
latter's position X2 being variable. Here xq = —0.01, m = 1, a = 
1/400 and xi = 0. 
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Figure 13: The left-hand Casimir forces in fig. 12 for a — 5 and X2 — 0.1 
plotted versus lattice constant a. 

In fig. 12 we plot Fj^ computed locally using {T^^){xi,X2,X3,a) and globally using 

. s _ E{xi,X2 + a,a) - E{xi,X2,a) , . 

Fl{xu X2) = . (4.30) 

a 

Again we find very good agreement as long as a is not too large and X2 not too 
small. For large a and small X2 the influence of V{x) is quite strong and we have 
to go to smaller lattice constant to observe the local-global agreement. We show 
in fig. 13 how this works for a = 5 and X2 — 0.1. Linear extrapolation to a = 
in fig. 13 yields the global and local values -56.073 and -56.026 respectively for Fl; 
again the agreement is extremely close. 

The final three figures in this section show Casimir forces computed globally using 
the vacuum energy. This is simpler than working with (T^^) which requires knowl- 
edge of eigenvalues and eigenvectors. However we have verified that (T^^) provides 
the same results. 

Fig. 14 shows how changing the position xq of the center of the Bessel boundary 
affects the Casimir force on boundary X2 in the system with three Dirichlet bound- 
aries at xi, X2, X3, all embedded in the potential V{x) — (a^ — l/4)/(a; — xqY. Fcas 
in fig. 14 is the net Casimir force on boundary X2 with xi, X3 fixed and variable xq. 
Three different X2-values are used. Essentially in fig. 14 we are holding all three 
Dirichlet boundaries fixed and rigidly translating V{x) by changing Xq. Clearly Fcas 
depends more sensitively on xq as a increases, which one would expect. Moreover 
the sensitivity of Fcas to Xq increases as the distance X2 — Xi between the leftmost 
and middle Dirichlet boundaries decreases. This is also to be expected. The con- 
clusion is, that one has to choose xq very small in order to obtain relevant results 
for Fcas- Our Choice xq = —0.01 seems to be appropriate. 

Finally we investigate the dependence on 0:3 of the Casimir force Fcas on boundary 
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Figure 14: The Casimir force ( |4.26| ) as a function of the position xq of the 
Bessel boundary for X2 = 0.5, 1, 2. Here a = 1/500, m = 0, a;i = 
and ^3 = 6. 
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Figure 15: Global Casimir force ( 4.26|) for variable position X3 of the rightmost 
Dirichlet boundary. Here xq = —0.01, m = 0, a = 1/500 and 
Xi = 0, a;2 = 1. The region a; > X3 is ignored. 
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Figure 16: Fr from eq. ( [4.29 ) for variable 0:3. Same parameters as fig. 15. 



X2- For large X3 this force (fig. 15) is practically independent of X3 as one would 
expect. At relatively small X3 a quite strong dependence on X3 abruptly appears 
and Fcas reverses direction - again quite predictable. The value of X3 at which Fcas 
vanishes (i.e. changes sign) moves inward with increasing a as it should. 



Fig. 16 displays the part of the Casimir force in fig. 15 coming from X2 < x < X3. 
For a > 1/2 this force component does not vanish as 0:3 ^ 00 (as it would if V{x) 
were zero). The Casimir force on the Dirichlet boundary at 0:3 obviously tends to 
zero as X3 — > cx) because V{x^) — >• 0. Thus the asymptotic force for large X3 in fig. 16 
is entirely due to V{x). Naturally this asymptotic force increases with increasing a. 



5 Conclusion 

Continuum QFT has always been used to formulate and study the diverse topics 
and systems of static and dynamical Casimir theory. We decided several years ago 
to reformulate a broad range of these topics and systems in the language of lattice 
QFT, confident this would provide a tool powerful enough to solve many Casimir 
problems not accessible to continuum methods. Here we have presented some re- 
sults on static finite lattice Casimir systems. A specific static spatial configuration 
was chosen: parallel Dirichlet boundaries with fixed separation on a d-dimensional 
lattice. We formulated this problem first with no background potential, then with a 
background lattice potential V{n) depending on the coordinate x = na perpendic- 
ular to the boundary planes at x = 0, /. Even for vanishing background potential 
V{n) = the presence of hard boundaries introduces technical features into the lat- 
tice formulation not present in lattice QFT with no boundaries. We have examined 
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many of these things in detail. With the introduction of a nonzero background po- 
tential V{n) there appear additional technical complexities which have also received 
careful attention. 

Renormalization of lattice QFT in the continuum limit a — is a prerequisite to 
success. One does lattice QFT primarily to be able to learn things about the contin- 
uum theory which are difficult to find out using continuum methods. Certainly that 
is our attitude toward this work. Until one has understood how renormalization 
works one's lattice results do not mean very much: the (informationless) a — > 
divergences will swamp any physical information contained in the lattice mathe- 
matics. For this reason we have paid close attention to renormalization, both for 
V{n) = and for V{n) > 0. Rather detailed results have been presented for one- 
dimensional lattices in sees. 3 and 4. We believe these results reveal a self-consistent 
and clear picture of the lattice QFT as it smoothly becomes (after renormalization) 
the continuum theory in the limit a — 0, both for vanishing and for nonvanishing 
background potential V{n). 

In unpublished work we have studied the extension of the V{n) = ID-lattice 
analysis to c? > 2 spatial dimensions. One can discern in the mathematics of these 
systems how lattice Casimir QFT in higher dimensions including renormalization 
parallels continuum theory. This is very important for us, because from the outset 
we planned to study 2D and 3D lattice configurations of hard boundaries which, 
because of nonsymmetric or otherwise complicated boundary geometry, cannot be 
dealt with using continuum mathematics. Work on such systems is in progress. 
Another direction of research we are pursuing is the introduction of background 
potentials V{n) on d > 1 spatial lattices. The semihard Bessel potential ( |4.2| ) is 
very useful in continuum theory because it leads to spatial modes which are known 
explicitly. However on the lattice it matters little whether one uses this potential or 
some other: The numerical work is made neither easier nor more difficult by one's 
choice of potential. Thus we are trying out on the lattice various semihard poten- 
tials to gain insights one would perhaps not be able to obtain differently. Another 
interesting type of enclosing boundary - soft boundaries - is also being studied on 
the lattice. We intend to report on different aspects of this work elsewhere. 
Beyond the static systems mentioned we have also been able to make progress in 
another important area: dynamical Casimir theory. One can, on a lattice, simulate 
moving boundaries quite effectively. Time-dependent backgrounds excite the quan- 
tum field and one can see this clearly in one's numerical results. We shall present 
lattice analysis of dynamical as well as static Casimir systems in subsequent articles. 
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A Operators on the lattice 



On a lattice one can reduce any differential operator to an ordinary numerical matrix 
once the fields 0„ have been labeled and ordered. Here we consider lattice versions 
of the following operators 

(j)A(j) andV{(f)V(j)) . (A.l) 

On a lattice there is, of course, more than one way to represent differentiation. We 
shall not be concerned here with distinguishing among different possibilities; the 
simplest and most direct definition will always be chosen. 



A.l Laplace operator 

The lattice Laplace operator is relatively straightforward due to its symmetry. How- 
ever, the precise form of the matrix representing A will depend on the boundary 
conditions imposed on the field (p. One can see this from the following two simple 
examples. 

Consider a ID lattice having five points labeled by i = 0, . . . , 4 and Dirichlet bound- 
ary conditions at i = 0, 4 ; thus the field consists of 



()^0,</^l,<)^2,</^3,()^4, 

with 



(A.2) 



00 = 04 = . (A.3) 
The second derivative is defined on the lattice in the usual way 

/ 9^0 \ 0i , 1 - 20, + 0i_i 



2 



(A.4) 



Because of the Dirichlet conditions we need - e.g. in eq. ( p. 411) - this operator only 
for the internal points riiii = 1, 2, 3). Thus 

3 

(A0), = ^Ai,,0j, z = l,2,3 with 
Ajj = — — 26ij + 5t-ij) , hj = 1,2,3 or in matrix form 



A = - 1-2 1 . (A.5) 




As a second example consider a ID lattice having four points and a periodic boundary 
condition; thus the field consists of 



^1, V2, (P3, (Pi 



(A.6) 
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The lattice point n = 5 is identified with the lattice point n = 1, i.e. 0i = 05 . With 
these boundary conditions we need - e.g. in eq. (|2.41|) - the second derivative also 
on the lattice "boundary", i.e. on the lattice point n = 1. One finds 



/-2 1 



1,2,3,4 with 



A = — 



V 



1 



1 



1\ 

2 1 

1-2 1 

1 -2 / 



(A.7) 



Using the periodic S function we can also write 



A 



1 



1,2,3,4, 



From these two simple examples one can immediately write down a general expres- 
sion for the Laplace operator appearing in eq. ( |2.41|) 



3=1 " 



(A.9) 



Here the ^-function 5„ „/ is appropriately periodic. This notation is, for several 
dimensions, neither very transparent nor very well suited for numerical work. One 
can however relabel the matrix elements according to e.g. the scheme 

(ni, . . . , nrf) ^ m, (A.IO) 
m = {ni-l){2N^Y-^ + {n2-l + N^){2N^Y-^ 

+ ■ ■ ■ + {ud-i -l+N^) {2N^) +nd + N^, 
1 < ^1 < iVi - 1 , -N^ + l<n,<N^ {i = 2,---, d) . 

This produces a symmetric matrix Am',m , whose dimension (A^^i — 1) (2A^^)'^~^ 
equals the number of lattice points. 



A.2 The operators (V^)^ and V(0V0) 

The lattice representation of the differential operator V is not unique. Let us again 
use the periodic ID lattice with four lattice points to illustrate this. We consider 
three definitions of V: 

(1) (V0). = ^^±i^ 



whose matrix representation (V0)i = Z]j=i -Dij 4>j is 



/ -1 1 \ 



D 



(1) 





1 



-1 1 
-1 1 
-1 / 



(A.ll) 



(A. 12) 
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(2) W)i 



(A.13) 



whose matrix form is 





f 


1 










1 




-1 


1 








a 







~1 


1 















-1 


1 y 



-{Si J - Si-ij) , 

(Jj 



(3) {V<f>)i 



2a 



(A.14) 



(A.15) 



whose matrix form is 



/ 1 -1 \ 

-10 10 

0-101 

VI 0-1 y 



(A.16) 



Definitions (1) and (2) have the disadvantage of producing nnsymmetric matrices. 
However, it is not our purpose here to discretize V0, but rather to define the lattice 
operator (V0)^. 

For periodic boundary conditions one has in the continuum theory 



J dx{V(t){x)f ^- j dx(t){x)A(t){x) 



(A.17) 



The analogous statement on the periodic lattice is (here k = 1,2,3 labels the three 
definitions) 



« E(E Di:j<t>i) (E - « E E 0KE(^^'^ ^knDl^M) ^ -a E 0,A,,>, . 

n i j i j n i,j 



It should be true that 

n 

However, when we calculate the left-hand side of this equation we find 



(A.18) 
(A.19) 



(2) 



(2-1 -1 \ 

-1 2-1 

0-1 2-1 

V -1 0-1 2 y 



(A.20) 



and 



(3) 
n,3 







2 





-2 


\ 


1 







2 





-2 






-2 





2 







I 





-2 





2 y 



(A.21) 
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Since we use the lattice Laplace operator ( |A.7| ) for periodic conditions we see that 



eq. ([A. 19|) is fulfilled only for the unsymmetric discretizations k = 1,2. Thus we 



are led to our final (symmetric) choice of (V0) 



1 



2 I \ a 



Wl-20„ + 0„-l , 0^+1 - 202 + 02_^ 



+ 



a2 2a2 

J,2 o^2 I J,2 



-0„ 2^ A„_„,(/)„/ H ^ — . (A.22J 



Using 



(V0)2 = V(0V0) - 0A0 (A.23) 
we obtain simultaneously from ( |A.22| ) a symmetrical definition for V(0V0) 

(V(0V0))„ = '^"-^^ ~ ^^^^ + '^-^ . (A.24) 

An important mathematical statement for periodic boundary conditions is the "spa- 
tial integral" 

/ dxV((f)V(p) Va(V(0V0))„ = 0. (A.25) 

^0 n=l 



With this lattice definition of /g dx eq. ( |A.17| ) is satisfied on the periodic lattice. 



Eqs. ( [A.22|) , ([A.24|) are also applicable for Dirichlet boundary conditions as long 
as one stays away from the boundaries. Boundary points have to be dealt with 
separately. For the left and right boundary points n = 0,N one can define right- 
and left-handed derivatives as follows 

(V0)2 = (V(0V0))o = 40?' (A.26) 
(V0)^ = (V(0V0))^ = 40^_i. (A.27) 



Then we can also use eqs. (|A.22 ), ( A.24|) at Dirichlet boundary points if we adopt 
the definitions 



0-1 = 01, (A.28) 
i^+i = 0^_i. (A.29) 



In the massive theory without potential we need more specifically </)_i 



(pN+l — —(pN-l- 

Including eqs. (|A.28|) , (|A.29|) we are led to a mathematical statement parallel to eq. 
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( [A. 251 ), but for Dirichlet boundaries 



/ N N-l 



2 



J2 a (V(0V0))„ + E « (V(0V0))„ = 

\n=l n=0 / 

N-l ^ ^ 

= « (V( 0V0))„ + — (V(0V0))^ + — (V(0V0))o 
= 0. 

Thus the appropriate lattice definition of "spatial integral" for Dirichlet boundaries 
is 

I 1 / ^ 

da: ^ -a{J2+Y.] , (A.31) 



/ 

Jo 



\n=l n=0 / 



and eq. ( |A.17| ) is again satisfied. 

In summary: For the basic lattice of this paper with Dirichlet boundaries at xi = 0,1 
and d — 1 periodic directions we have, with eq. 



ml = -<^>n E ^n,n'^n' + E "^^^ Z; ^""^- , (A.32) 



n/6G j=l 



2a2 



(V(0V0))„ = E ' o 2 ^"^^ • (A-33) 



B A lattice artefact 

In calculating the {d = l)-vacuum stress tensor in sec. 3 with V{x) = we encoun- 
tered a lattice artefact in the lattice formulation which (i) diverges in the continuum 
limit a — i> and (ii) has no counterpart in the continuum theory. The origin of such 
a lattice quantity is to be found in the noncommutativity of the mode sum over 
"momentum" k with the limit a — 0. Our task in this appendix will be to establish 
some needed properties of this lattice artefact term 

/v-i 2iv2 sm\^) cos(^) , / , , 

A(mJ,n,a) = -Y V2iv^ v n i ^^^^ ^ ^_ 

fc=i P^{mlY + m^ sin^dl) a 

Most importantly we will show that A contains no physical information. 
To illustrate noncommutativity let us observe that 

lim A(m,/,n,a) 7^ (B.2) 

a— >0 



51 



while in the reverse order 



Note that A has no global counterpart 



N 

^ A(m,/,n,a) = (B.4) 

n=l 



because J2n=i cos{27rkn/N) = for /c 7^ 0. Expanding A in powers of m we find 

^ , , , 1 . vrA; 27rkn. 

A{mJ,n,a) = 2^ sm (^) cos(^^) 

m^^M . .7Tk 2Tikn 
H > — sm cos B.5 

8 ^2N' ^ N ' ^ ' 

'"^-1 cosr^-'^"^ 



N 



l'^ u=2 V 2 / k=i m (A^sin(^ 

I + II + III . 



2N 



2u-3 



The term III : 

One can verify (numerically for example) that for Re > —1 

/^-i cosf^^) \ 
lim y V jv ; 1 Q Vn,Re/i>-l. (B.6) 

In the contribution III we need this statement for /i = 2z/ — 3 > 1. To illustrate the 
limit (BJ) consider the /U = 1 case. Because 

2 TT 

sin(x) > —X, < X < — , (B.7) 

TT 2 

it follows that 

V '"'^ ^ ^ <^v|<- = 0(a) Vn. (B.8) 
iV2iVsin(i|)-iV2^^^fc-iV ^ ^ 

For yU > 1 

because the sum over k now converges. In summary, the contribution III is at most 
0{a) and therefore this contribution to A vanishes in the continuum limit. 
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The term II : 
We define 

^tzl 1 / vrfc , ,2'nkn. , ,1 ^ ,^ ^ 

fc=i 

Clearly 

^ I 

J2hi-,n) = 0. (B.ll) 

n=l 

The finite sum ( |B.10| ) can be evaluated in closed form 

^(a'"") " 4/ (tan(f + vrf ) + tan(f - vrf ) " ^) ' ^^'^^^ 
Writing n = x/a and using the formula 

^ + -^^ = -2(1 + -^^) a + 0{a^) fiir a < x (B.13) 



tan(a + x) tan(a — x) V tan^x 

we obtain 



h (L -]=-^ + 0{a^) , ^ « X « (/ - ^) . (B.14) 
\ a a I 21 4 4 



In the continuum limit 

lim/i ( -, - I = 0, 0<x</ (B.15) 

and therefore 

, A ^ , / 1 , ttA; , ,2T:kn\ , 

lim/i -,- = y lim — sin — cos — — = 0, 0<x</. B.16 

For < X < / it does not matter how one takes the limit a — > 0. However on the 
boundaries x = 0, / we have 

\a J \a J 21 Vtan(f ) / vr ^ ^ 

In other words, the (A^ oo)-limit of the term II is nonuniform in lattice position 
n. This phenomenon is quite familiar in infinite series. We see in eq. ( p.l2|) that h 
has poles at n = 1/4 and n = N — 1/4. These poles lie between the boundary and 
the first adjacent lattice points - i.e. they are not on the lattice. They signal to us 
that the lattice function h is nonuniform as the boundaries are reached. 
According to eq. ( |B.16 ) the continuum limit of II is zero between the Dirichlet 
endpoints. At these endpoints the value of II in the continuum limit jumps to 
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m?/4:Ti. As there is no / dependence in this contribution to A one sees that the 
Casimir force is unaffected. 
The term I : 
We define 

fc=l 

where the function / can also be written 



The finite series (B.18) can be evaluated in closed form 



1 ■i.Tck. ,2'Kkn. 1 ,^ , 

— y sin^ cos = — B.20) 

al ^2N' ^ N ' 2al ^ ^ 

(I 1 3 2> \ 

+ 



16a/ Vtan(|| + 7rf) tan(|| - vrf ) tan(^ + vrf ) tan(3^ - vrf ) 
and 

, . _ 1 / 1 1 



16a/ V tan(^ + tt^) tan(^ - ^T'-^ 
3 3 



tan(^ + 7r— ) tan(^ - vr— J 



(B.21^ 



We see from eq. ( [B.21| ) that /(/, n, a) has poles atr2 = l/4, 3/4 and n = / — 1/4, / — 



3/4. These poles lie between the Dirichlet boundaries and their first adjacent lattice 
points. From this we anticipate that the function / will also exhibit nonuniform 
behavior as the boundaries are reached. 



From eq. (|B.13|) it follows that 



f{l,-,a) = 0{a'), ^«x«/-^. (B.22) 
a 4 4 

This means that / vanishes in the limit a —>■ everywhere between the boundaries 

lim/(/,^,a) = 0, 0<a;</, (B.23) 



and that for < a; < / we can in eq. (p.l9| ) freely interchange J2k with (a —>■ 0). 



The behavior of / on the boundary is obtained from the expansion (for n <^ N) 

fa n,a) = ^( + — ^) +0(a') , n <^ N . (B.24) 

' ^ 4a2V3 + 4n 3-4n 1 + 4n 1 - 4ny ^ ^ ' ^ ' 

In the limit a — this expression diverges; however the value approached is in- 
dependent of /. The Casimir force on the boundary is the discontinuity in (T^^) 
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across the boundary. An /-independent jump in (T^^) makes no contribution to the 
Casimir force. On the boundary J2k s-nd a — > are not interchangeable. Note also 
that 

N ^ 

Y^af{l,n,a) = --. (B.25) 

n=l 

In summary we have for the lattice artefact the following: 

1 fn?' 

A{m,l,n,a) = — + f{l,n,a) + —h{l,n,a) + 0{a) , (B.26) 
zcil 8 

N 

Y,A{m,l,n,a) = 0. (B.27) 

n=l 

A{m, n, / , a) — 1/ (2a/) vanishes in the continuum limit a — >■ for every point between 
the Dirichlet boundaries. This function on the boundaries x = 0,1 jumps to a 
nonzero value which - being independent of / - cannot influence the Casimir force. 
The comment now to follow is speculative and mathematical. Let us define a lattice 
function 

6(1, n, a) := -2./(i, n,a) = U - ' . " wS^T^) <^-''' 



tan(- 

3 



tan(f + 7rf) tan(f - vr^f ) 



which seems to be rather like a 6 function. Eq. ( |B.25|) now reads 



N 

^a5(/,n,a) = 1. (B.29) 

n=l 

The properties of / discussed above show that 

, . / 0(a) <n<l m Qn^ 

= I 0(l/a) n = 0,l (^-30) 

and the "integration" (|B.27D over the lattice only receives a contribution from the 
boundaries. This suggests that 

lim 5(/, n,a) = l {6{x) + 6{x - /)) . (B.31) 

a— >0 2 

To test this numerically we applied a test-function polynomial t{x) = Ax^ + Bx^ + 
Cx + D and found that, with very high accuracy 

1 / ^ iV-l \ ^ 

- ^a5(/,n,a)t(na) + ^ a(5(/,n,a)t(na) ^ - (t(0) + t(/)) , (B.32) 

2 Vn=l n=0 / 2 

712 (a) 

aS{l,n,a)t{na) ^ 0. 

n=ni (a) 
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In the latter sum, for each lattice constant a used we required that 

ni{a) a = const = Xi > and Xi < ^2(0) a = const = X2 < I ■ (B.33) 
Finally we mention that 

21 \tan{^ + vr— ) tan(^ - TC—jJ 
1/1 1 \ 



(52(/,n,a) _ 2/ (^tan(3f + 7rf ) +tan(3f -7r^)y 

have the same properties as 6{l,n,a). Also 

3 1 

S{1, n, a) = a) + ^'^2(^, a). (B.35) 

C Perturbation theory on the one-dimensional lat- 
tice 

Here we calculate to first order in the lattice potential U{n) = + V{n) the lat- 
tice vacuum energy E^ac and stress tensor component (T^^) for the one-dimensional 
lattice with Dirichlet conditions at its endpoints. The "unperturbed" system is the 
massless scalar field with vanishing background potential V{n) = 0. Our calculation 
for arbitrary U{n) will determine the additional divergences in E^ac and (T^^) arising 
from U{n). 

For a given lattice potential one begins by first determining numerically the eigen- 
vectors and spectrum of the matrix 

On,n'{X) = An,n' + \U{n)Sn,n' (C.l) 



a2 



+ XU{n)6n,n', I <n,n' < N -1 



Then from these ingredients one can find the first-order shift in the vacuum energy 
and in (T^^). 

For A = the spectrum ul and eigenvectors w'' of ( |C.1D are known 

- (C.2, 

where 1 < < — 1 and / = Na. These eigenvalues are not degenerate, so we can 
use nondegenarate perturbation theory. 



56 



For A 7^ we name the eigenvalues and eigenvectors of ( |U.1| ) £:|(A) and {/^''(A) re- 
spectively 

N-l 
n'=l 

Expansion in powers of A yields (see e.g. ref. [|TB[) 

4(A) = u;l + Xel^'^ + X^el^'^ + ... , (C.5) 
t;'=(A) = ^i;'= + At;^W + AV'=(2) + --- (C.6) 



where 



-2(1) 



iw'Uw')=j:^^sin\^)U{n) (C.7) 



/ sin^( — — )f/(x)(ia;. 
Jo I 



V 



I Jo 

p/fc 



Ground state energy: 

In the expansion of the lattice vacuum energy 

N-l 



iv — i -| 

E{X) = E Ae,(A) =E(°) + AE« + A2e(2) + --- . (C.9) 
fc=i 2 

the leading unperturbed term is 

To find the first correction to this we need the connection between the terms in 

ekiX) = M{X) =u;u + X si'^ + X' + ■ ■ ■ (C.ll) 



and the terms in eq. ( |C5| ). From the identity 



OO \ 2 OO 



E«-^' = Jlbux" (C.12) 

\iy=0 ) u=0 



OO OO OO OO I V 

u=0 u=0 ^l=0 u=0 \Ai=0 



one finds 



au = ^K^^h^,. (C.13) 

n=o 
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Thus 



z/ = : 
u = l: 

z/ = 2 : 



bo 
bi 



ao ■ 
ai 



(C.14) 



1 



a2 



4ao, 



Then from 



.(1) 



we find 



N-l 



2uk 



N-l 



(C.15) 



^ \p i (1) = 1 V 



sm 



fc=i 



,7Tkn 



) t^H • 



(C.16) 



This is the first-order correction to the lattice vacuum energy for arbitrary lattice 
potential V{n). 
Let us define 



with 



U{x) = U + U{x) 



U = \ [^U{x)dx. 
I Jo 

U{x) dx = . 



(C.17) 



(C.18) 





Substitution in eq. (|C.16|) yields with the help of (|3.15| ) 

N-l 

U [x) dx 

lo 



If 

8 io 
1 



'2 



J N-l 

E 



E^-^(^)^(-) 



I U{x) dx f- In + c{N) \ + - V 

Jo Vtt / 8 rr; 



dh 



where 



k=i Nsiny^^ 



- f sin^f — —)U{x)dx. 



(C.19) 



n=l 



/ Jo 



I 



(C.20) 



lim dk = . 

k^oo 
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The second term in the second equahty of eq. ( |C19| ) is finite for a ^ 0, while the 
first term is logarithmically divergent. 

Now we have what we need write down the divergent terms in E{X) 



21 

na 

l2 



_ _ ^ _ A 

2 ~ 2^ ~ 4^ 



In a / dxU{x) + terms finite for a ^ . (C.21) 
Jo 



Terms of 0{X ) and higher do not contain a — > divergences. This can be checked 
numerically. For U{x) = m? we recover eq. (|3.17|) . Note that on the lattice one can 
replace 



U (x) dx 



1 



/N-l N 



(C.22) 



\n=0 n=l/ 

The difference is at most of 0{a) and, because a Ina — > in the continuum limit, 
this difference is irrelevant. 

For the lattice (T^^) we have from eq. ( |2.77|) 



fc=i 



2aek{\) 



k=l 



From eq. ( |3.33| ) the unperturbed tensor is 



(C.23) 



(C.24) 



where / is the function ( |B.19D . The first order correction to this is found by substi- 
tuting eqs. (|C;6D and (^[11]) into 



^11 \ (1) 



N-l 




E 

fe=i 


2a 


N-l / 


+ E - 




k=l \ 


+ 





Si 



.k \2 



£k \ -\ W„'V,^ 

a 

k \2 



(C.25) 



(<-+lr-2K0^ + K.„l)^(l) 



Next we introduce the potential 

U{n) = U{n) - U{n') , i.e. U{n') = 



(C.26) 
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which leads to 



el^'^ = el^'^ + U{n'), (C.27) 



.(1) _ Ai),Uin') 

k — fcfr "T 



k{i) ~k{i) 



Here the quantities with tildes are defined with V in place of V . Continuing with 
this tilde notation we find for the first-order correction to (T^^)^'^-' 

(1) =ri \ U{n') 1 ^ U{n') ^1 1 TTfc 27rfcn' 

where h{N,n') is the function eq. (p.lO| ). 

The lattice quantity (Tl^^)^^^ remains finite in the limit a — and is the first-order 
correction to (T^^)^^^ for a lattice potential V{n) which vanishes at the point n = n'. 
It is therefore relatively obvious that {T^})^^^ contains no a — > divergent terms 
because the lattice should be (in first order) proportional to V{n'). For V{x) = 
eq. ( |C.28|) reduces to the system discussed in sec. 4.2. From (|U.24| ), ( |C.28| ) we can 
read off the divergent terms in eq. ( |C.23| ) 



(T„ii)(A) = ^ + /(/,n,a) + ^lna + ---. (C.29) 



Second-order corrections need not be discussed. Numerical probes show that all 
divergent terms are displayed in eq. ( P.29|) . 
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